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ABSTRACT 


The contribution of this thesis is the development of a CAD (computer aided 
design) tool for current mode multiple-valued logic (MVL) CMOS circuits. It is 
only the second known MVL CAD tool and the first CAD tool for MVL CMOS. 

The tool accepts a specification of the function to be realized by the user, 
produces a minimal or near-minimal realization (if such a realization is possible), 
and produces a layout of a programmable logic array (PLA) integrated circuit that 
realizes the given icaear The layout is in MAGIC format, suitable for submission 
to a chip manufacturer. ‘The CAD tool also allows the user to simulate the realized 
function so that he/she can verify correctness of design. 

The CAD tool is designed also to be an analysis tool for heuristic minimization 
algorithms. As part of this thesis, a random function generator and statistics gath- 
ering package were developed. In the present tool, two heuristics are provided and 
the user can choose one or both. In the latter case, the better realization is output 
to the user. The CAD tool is designed to be flexible, so that future improvements 
can be made in the heuristic algorithms, as well as the layout generator. Thus, 
the tool can be used toaccommodatenew technologies, for example, a voltage mode 


CMOS PLA rather than the current mode CMOS currently implemented. 
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I. INTRODUCTION 


A computer-aided design (CAD) tool has been created at the Naval Postgrad- 
uate School to design multiple-valued logic programmable logic arrays (MVL-PLA). 
The CAD tool accepts a function specification and produces a circuit layout. Small 
MVL-PLA’s can be designed by hand. However, for most practical applications, 
hand design is too time consuming and prone to errors. 


The CAD tool performs the following: 


e Accepts an input function specification, specifically an algebraic expression. 


e Minimizes the expression producing a realization that occupies the least 


amount of space. 


e Produces a layout which can be used to generate the MVL-PLA chip. 
The CAD tool is designed toaccommodate future improvements in: 


e The minimization program and 


e MVL-PLA layout. 


It is expected that as knowledge is gained about these aspects of MVL-PLA 
design, there will be new subroutines that will substitute for these parts of the tool. 
The whole program has been written so that such changes can be easily made. 

This CAD tool, which applies to current-mode CMOS technology, is only the 
second known tool for multiple-valued logic circuit design. A CAD tool exists for 


the design of multiple-valued CCD PLA design [Ref. 1]. 


A. BACKGROUND 

The promise of current-mode CMOS circuits for compact realizations has gen- 
erated much interest. Japanese researchers have developed a binary multiplier in- 
tegrated circuit using bi-directional current-mode multiple-valued CMOS [Ref. 2]. 
The circuit has speed almost equal to the fastest binary multiplier [Ref. 3] but 
requires half the chip area and dissipates half the power. Researchers in Holland, 
France and the United States have also fabricated current mode multiple-valued 
CMOS circuits. 

Since this technology is new, no design tool has been developed for it. Design 
tools are absolutely needed to produce practical circuits since modern VLSI circuits 
are so complicated. It is the lack of such tools that has motivated the research 


described in this thesis. 


II. PROGRAM DEVELOPMENT 


The package developed in this thesis is a comprehensive program written in 
the “C” language to run on the VAX-11/785 computer at the Naval Postgraduate 
School. As “C” is highly structured, the program is easy to maintain and easy to 
develop into a large program package. Also, use of “C” assures compatability with 
presently available binary VLSI design tools such as the Berkeley CAD tools [Ref. 
4]. 


A. APPROACH TO THE PROBLEM 


A top-down programming technique was used. First, the tasks were defined 
which the program needed to perform. Then, each task was developed through a 
series of steps to identify very specific routines required to perform each task. 


The major tasks are as follows: 
e Data Entry 
e Logic Minimization 


e Logic Simulation 


Layout Generation 


1. Data Entry 


Allows the user to enter input expression from an input data file. 
— Partitions input expression into identifiable parts, called tokens. 
~ Examines each token for errors. 


— Flags errors when they occur. 


2. Logic Minimization 
— Three minimization algorithms are used. 


* Pomper and Armstrong [Ref. 5]. 
* Dueck and Miller [Ref. 6]. 


* Gold [Ref. 7]. 


All of them are heuristic and based on the direct cover approach. 


Minimization produces a minimal or near miniinal solution, thus reducing 


the chip area needed. 


This routine will be discussed in more detail in the following section. 


3. Logic Simulation 


— Verifies the minimal solution from the Logic Minimization routine by 


applying test vectors given by the user to it. 


4. Layout Generation 


— Uses a minimal solution from the Logic Minimization routine to produce 


a layout which will be used for MVL-PLA chips. 


B. ORGANIZATION OF THE PROGRAM 

The basic organization of the program is shown in Figure 2.1. The solid lines 
show flow of control. The user supplies input expressions via the data entry routine. 
These are applied to the logic minimization routine for obtaining the minimized 
expression. The minimized expression is then passed to the logic simulation routine 


for verification. Then it is used for layout generation for MVL-PLA chips. 


A CAD TOOL 


Layout for 
Input MVL-PLA 
Expression 
Logic igen Chips 





Minimization Simulation 





Figure 2.1: Program Organization 


PLA design using this program resembles program development under an op- 
erating system. An editor is used to develop the source file (input expression). ‘This 
is compiled by the data entry routine. If it 1s free of syntax errors, it is assembled 


(minimized) and run (a layout is produced). 


Ill. PROGRAM DESCRIPTION 


The CAD tool is best described by considering a specific example. The user 
starts with a logic diagram (or Karnaugh map) of the input expression to be realized. 
An example of the logic diagram is shown in Figure 3.1. It shows a 4-valued 2 
variable expression. Without a CAD program, an inexperienced user might spend 
half an hour to find a minimal solution of the expression. What is worse, when 
the expression is complicated because of a large radix and/or a large number of 
variables, it may be ianeeeiule not only to find the minimal solution, but also to 
visualize it. 

Since the CAD program ts designed to find a minimal solution of the expres- 
sion, the user is not required to mimumize the expression. All the user has to do is 
specify the expression correctly. Figure 3.2 shows one way the function of Figure 
3.1 can be expressed for the CAD tool. Here each circle represents a product terin. 


Equation 3.1 shows an algebraic equivalent of the function of Figure 3.2. 


Ars Git) 
+2 +* X1(1,2) + X2(0,0) 
“3 +X 1333) exe gie 
“2% X ICO Oe 2 (eZ) 
+1 * X1(0,0) + X2(2,3) 
+2 + X1(1,1)+* X2(1,1) 


+2 * X1(1,1) * X2(2, 3) 


+1 * X1(1,1) * X2(2, 2) 
42 + X1(2,3) * X2(1, 1) 
ele 2 2X 2(1, 2) 
+1 * X1(3,3) * X2(2, 2) 
+3 *.V1(2,2) * X2(3, 3) 
Equation 3.1 1s represented according to the standard format of the expression 


shown in Table 3.1. Any multiple-valued logic function f(X1,.X2,X3,---,Xm) can 


be represented by thie followmg format: 


i radix : number of variables : ee), 
+c, * X 1(al, 61) * X2(a2,b2) *--- * Xn(an, bn) 


+¢e, * Al(cl,dl) * X2(c2,d2) *--. + Xn(en, dn) 


+c, * X1(f1,g1) * X2(f2, 92) *---* Xn(fn, gn) 


where c},C€2,::: and c, denote a logic value between ‘1’ and ‘(radix -1)’, the symbol 
‘*’ denotes the minimum operator, Xn(an, bn) denotes a literal function of the input 
variable ‘Xn’ with lower bound ‘an’ and upper bound ‘bn’. Specifically, Xn(an, bn) 
is (radix -1) if an < Xn < bn, and 0 otherwise. The symbol ‘+’ represents the 
truncated arithmatic sum (truncated to (radix -1) if the actual sum exceeds (radix 
a) 

According to Equation 3.1, the expression has radix 4 and 2 variables and a 
list of terms (11 terms in this expression). A term is composed of a coefficient and 


a list of variables (2 variables for each term in this case). The expression shown 11 


Equation 3.1 is entered into the CAD system via the data file. 


1 


xl 


NO 
Ce 
CO 
= 
it 


Figure 3.1: An Example of Logic Diagram for 4-Valued 2 Variable Ex- 
pression 





Figure 3.2: An Example of Logic Diagram for An Input Expression Spec- 
ified by the User 


Expression Radix : Number of Variables : A List of Terms 


Coolliclent * A List of Variables 


Varlable X ( Vartable ID ) ( Lower Bound, Upper Bound ) 





TABLE 3.1: Standard Format of an Input Expression for a CAD Tool 


A. DATA ENTRY 

Figure 3.3 slows the flowchart of the Data Entry routine. The input expression 
from the input data file is applied to this routine first. It examines all the expressions 
in the data file. There is a special function, called the parser/scanner in this routine. 
It partitions the expression into identifiable parts, called tokens. Then, it passes each 
token to the next three consecutive error checking subroutines — the syntax error 
checker, the grammar error checker, and the range error checker. 

If an error is found, the corresponding error message is presented to the user 
and the program execution stops. The user needs to go back to the input data 
file and correct it. If all the expressions are correct, all are stored in the allocated 
memory space. Figure 3.4 shows the data structure of the input expression. It is 


developed to accommodatethe input expression. 
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Figure 3.3: Flow Chart of the Data Entry Routine 
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Figure 3.4: Data Structure of the Input Expression 


IZ 


According to the Figure 3.4, three data structures form a complete input 
expression. They are the Expression, Term and Bounds data structures. The radix, 
number of variables and number of terms in the input expression are saved in the 
Expression structure. 

The size of the array for the term lists is determined by the number of terms of 
the expression. Also, the size of the bounds array is equal to the number of variables 
used in the expression. Therefore, all the components of the input expression are 
saved into this data structure. 

Once this data structure for the input expression is formed, it will be used 


throughout the CAD program. 


B. LOGIC MINIMIZATION 

The data structure of the input expression is formed in the Data Entry routine. 
This is then used by Logic Minimization routine to minimize the expression. Figure 
3.5 shows the flowchart of the Logic Minimization routine. 

First, a working copy (f) of the original input expression is made. If f is fully 
covered (made up of only zeroes and don’t care terms (radix)), the program stops. 
The minimization process is accomplished. Otherwise, the prograin finds the most 
isolated minterm first, then picks the best implicant (or product term) from among 
the possible implicants which cover the chosen most isolated minterm. 

The best implicant chosen is saved as a part of the minimal solution of the 
input expression. Then it is taken out from the working copy (f). This process is 
repeated until the resulting expression has no more minterms. 

As mentioned earlier, three heuristic algorithms are used in this thesis. As 
shown in Reference 7, neither that of the Dueck and Miller [Ref. 6] nor that of and 


Pomper and Armstrong [Ref. 5] is consistently better than the other. That is, there 
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Figure 3.5: Flowchart of Logic Minimization Routine 


are functions where the first does better than the second and vice versa. Thus, the 
technique of Gold is introduced in Reference 7. This is simply the application of 
both algorithms followed by a choice of the best realization. Since the CAD tool is 
also used to analyze minimization algorithms, implementation of all heuristics is an 
integral part of the program. All of them work in an identical manner except in two 
ways. One is the way a minterm is cllosen, and the other is the way an implicant is 


chosen. 
e Pomper and Armstrong [Ref. 5] 


The minterm is chosen randomly. However, tle implicant is cliosen as the 
one which drives the most minterms to zero when subtracted. I{ more than one 
implicant is available, the largest is chosen. If there are more than one largest, the 


first generated is chosen. 
e Dueck and Miller [Ref. 6] 


In this algorithm, an isolation factor (IF) is calculated for each minterm with 
the smallest value beginning with all 1 minterms, 2 minterms, etc. The minterm 
with the largest IF is chosen and if there are more than one, the first generated is 
chosen. IF is calculated as the inverse of the clustering factor where the clustering 


factor is as follows: 


Ch — DWiA. x \nages — Weer Ae (GE). 


where DEA is the number of variables (directions) in which a minterm (a) can 
be combined with a non-zero number of minterms (called direction of expandable 
adjacency). EA is the number of adjacent minterms with which a minterm (a) can 


be combined in an interval literal (called expandable adjacency). CF provides a 
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measure of the degree of which a specific minterm can combine with other minterms 
in the expression. 

The minterm with the smallest CF is chosen as the most isolated minterm. 
All implicants that cover this minterm are then generated. A measure called the 
relative break count (RBC) is calculated for each. This provides a measure of the 
degree to which the expression is simplified if the implicant under consideration is 
chosen. The implicant with the smallest RBC is chosen as the best implicant. If 


there are more than one best implicant, the first generated is chosen. 
e Gold [Ref. 7] 


Gold is a heuristic in which the algorithms of Pomper and Armstrong and 
Dueck and Miller are apphed, and the best realization (the one with the smallest 
number of product terms) is chosen. It was inspired by the observation that these 
algorillims display a diversity in realizations and one algorithm is not consistently 
better than the others over all expressions. However, certain classes of expressions 
do show that one algorithm does perform consistently better. With Gold, the com- 
bination of algorithms takes advantage of the best features of each. 

Table 3.2 shows the comparison of three heuristic algorithms. When logic 
minimization is accomplished, the minimal solution of the original input expression 
is as shown in Figure 3.6 and Equation 3.4. Note that the number of product terms 
in the input expression in Equation 3.1 was 11 (eleven). Now the minimal solution 
(Equation 3.4) has 8 (eight) product terms in it. There is about 25% reduction in 


the number of product terms. 
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Figure 3.6: The Minimal Solution of Input Expression from Logic Mini- 
mization Routine 


Choice of | Chovecwen 





Alyorithin ) 
Minterim Implicant 
4. Pomper and Drives Most 
Random 
Armstrong Minterms to O 
Largest IF Smallest ABC 
2. Dueck and 
(isolation (relative break 
Miller factor ) count ) 
3. Gold Best of 1 and 2 


TABLE 3.2: Comparison of Minimization Algorithms 


ie (3.4) 
+1 * X1(0,2) + X2(2,3) 
+] 4X1 (3, 3) ox ie) 
+1 * X1(1,1) *.X2(3, 3) 
+] X12, 3) eee le) 
42 * X1(2,2) * X2(3, 3) 
+2 * X1(1,3) * .X2(0,0) 
2 xX (ORZ ees 2 ll) 


42 + X1(0,1) * X2(2, 2) 
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C. LOGIC SIMULATION 

This routine is used to verify the minimal solution from the logic minimization 
routine. The user can specify certain specific values or range of values for each 
variable in the input expression (called test vector). The logic simulation routine 
evaluates the minimized function value for each coordinate from the test vector. 
The user can do the verification test by comparing the original expression with the 


output of the Logic Simulation routine. 


D. LAYOUT GENERATION 


The output of the Logic Minimization routine is used to produce a layout file 
(called a magic file) for MVL-PLA chips. This routine is discussed in more detail 


i Refereice 8. 


iL, 


IV. CONCLUSION 


This thesis presents a CAD tool for design of multiple-valued current-inode 
programmable logic arrays. It was designed to be user friendly, requiring little 
knowledge of programming techniques. Appendix A contains a user manual written 
for the first time user. In the interest of clarity, details of the program have been 
onitted. However, for the interested user, the program las been extensively docu- 
mented. This documentation and the program are contained in Appendix B. Three 


main conclusions have resulted from this work. 
e Conclusion 1 - Design/Analysis Tool 


The CAD tool can be used for two purposes. One ts as a design tool and 
the other as an analysis tool. It is used for designing practical current-mode MVL 
CMOS circuits and also used for analyzing and comparing the different minimization 
algorithms. In fact, it is designed to get the statistics for each minimization program, 


which is very useful to obtain insight into logic minimization algorithms. 
e Conclusion 2 - Further Improvement in CAD 


Further research and improvements are needed in two specific areas. ‘These 
are 1) Logic Minimization and 2) Layout Generation. Still there are cases which 
the heuristics give poorer realizations than supplied the user. Programs are desired 
which can find the minimal solution closest to the optimal solution. In MVL-PLA 
layout, there is some waste of the chip area. The further optimized layout configu- 


ration should be considered. 
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e Conclusion 3 - Execution Speed Improvement 


The Profiler Utility shows how many times each routine is called and the 
percentage of time spent in executing that routine. Below is an example output from 
this routine. According to the profiler output shown in Table 4.1, there are several 
subroutines which are called most frequently throughout the program execution such 
as .eval, _vcopy , and _next_coord. Typical runs show these three routines as the 
most time consuming. If those subroutines are carefully examined and improved, 


the total execution time will be reduced significantly. 


A) 


time cumsecs 


141 
143 
147 


wo 
geaoeoooooocoo0ocoo0o oa 0 co ooo aogag oa ONFN 


* 


147 


147 


147 


147 


147 
147 


147 


OSOoo0o0o0ocoa9#oecdoc‘coodoooqoqcoocoqo0ooqoao od oO B&B WO 


{6 
RES, 
-08 
147. 
147. 
1474 
147. 
147. 
147. 
147. 
147. 
20 
147. 
147. 
220 
.20 
147. 
220 
147. 
147. 
147. 
a2 
.20 
i Be es 
147. 
.20 
Lat: 
147. 
147. 
147. 


eS 
16 
17 
18 
18 
ie, 
20 
20 


20 
20 


20 


20 
20 
20 


20 
20 


20 
20 
20 
20 


fcall 
863184 
863184 
443878 
6133 
1789 
391 
358 
358 
S56 
244 
244 
TA0 
164 
162 
162 
158 
140 
iZ 
96 

80 

80 

80 

80 

pe) 

38 

32 

32 

a2 

a2 

32 


ms/call 
.16 
-00 
~O1 
.O00 
-O1 
.03 
hs) 
.00 
.03 
.04 
-00 
.00 
.00 
200 
.00 
00 
.00 
.00 
.00 
.00 
~J00 
ONG, 
200 
E00 
200 
ZOO 
200 
.00 
200 
0.00 


QO OOOO GOO O00 O00 OOo OO C0 0 OC So Oo oara Ge 


fhame 

me Gel 

_vcopy 
_nexticoond 
_next_ implicant 
_copy implicant 
_yylook 

_yylex 

next Coken 
TSE Ecae 
_streomp 

_match 

Jstrien 

fi 1puk 
__doscan 
_sscanf 
_malloc 

2 ee 
_alloc_bound 
_limit 

_alloc texm 
_var_id 

_var list 
_variable 
_realloc 
_becopy 

_gen bounds 
_init implicant 
_subtract_ implicant 
Scerm 

te Empelansie 


TABLE 4.1: The Output of the Profiler 


yap 


APPENDIX A 
Getting Started 


A few examples are included in this appendix A in order to demonstrate 
the use of this CAD tool,the appearance of a terminal session, and the 
capability of the utility functions. 


1. DESIGN PROCEDURES 
(a) LOGGING IN 
- LOG IN ON A VAX TERMINAL. 
- CHANGE TO SUBDIRECTORY 'mag' IN YOUR WORKING DIRECTORY. 
% cd mag 
(b) CREATING THE INPUT DATA FILE 
- BY USING ‘vi'EDITOR, MAKE AN INPUT DATA FILE. 


% vi tt2.mvl 


(1) WHEN THE USER HAS A SINGLE INPUT EXPRSSION 
4:2: 

ea XCl, 3) 7 X2(1719 

tar xl, 3) *X2(1,3) 

+2*X1(1,2)*2200,2) 

Ex Cl, 2)* X23, 3) 

+2*X1(0,2)*X2(0, 2) 

Tole, 3) * X2(2,3) 

#3*X1(0,2)*X2(0,0) 

date Xe 1, 2)*X2(0,2) 

+41*X1(0,3)*X2(1,2) 

tut X1(172)* X2(1, 2) 

+1*X1(1,2)*X2(0,1); 

( tt2.mvl is the name of input data file. In this example there is 
only one expression. It is a 4-valued 2-variable expression with 
11 terms in it. A ;(semicolon) terminates the expression. ) 


(2) WHEN THE USER HAS MULTIPLE EXPRESSIONS, THE INPUT FILE MAY LOOK AS 
FOLLOWS. THERE ARE THREE EXPRESSIONS AND ALL OF THEM ARE 4-VALUED 
2-VARIABLE LOGIC EXPRESSIONS. 


4:2: 
+2*X1(3,3)*X2(1,3) 
OP XTC1,3)*X2(1,2) 
Ve Xi(1],3)4%2(1,1) 
+1*X1(1,3)*X2(1,3) 
Boel, 2)*X2(0,2)1 
Yea 

tue ee 1, 2)*X2(3,3) 
t2*XU(0,2)*X2(0,2) 
pees ,3)*X2(2,3)7 » 
4:2: 

pie XT 1, 2)*X201,2) 
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(3) 


(2) 


(Cc) 


+E XE CZ) tk 2a 
+27 X00, 0) 320) 
+3*%1(0,1)*xX2(1, 3) 
(input data file; tt3.mvl ) 


WHEN THE USER APPLIES TEST VECTORS TO EACH VARIABLE FOR VERIFICATION, 
THE INPUT DATA FILE LOOKS AS FOLLOWS. 


Ac 

P2eR CI 3) exo Ie) 
+1*X1(1,3)*X2(1,3) 
#94X0(1, 2) oko) 
+1*X%1(1,094%209,4) 
4+2+*X1(0,2)*X2(0, 2) 
+ 2° X103 SO A227) 
Foe NCO 2) AZ Oe) 
ti xi(1,2)402 00,2) 
#1*X1(0,3)*X2(1,2) 
#1*X1(1,2)*X2(1, 2) 
+1*X1(1,2)*X2(0,1); 


£109,353) %200,.3); 
(test vector for the input expression above. Be sure to put ';' at the 
end. No need to put '*' between variables. In this case test vector 
has the range 0 through 3 for both variables X1 and X2. The user can 
specify any arbitrary test vector to each variable, for instance, 
X10, D202 2 ye 
( input data file; tt4.mvl ) 


WHEN THERE ARE DON'T CARE TERMS IN THE EXFRESSION,USE '?' (QUESTION MARK) 
FOR THE COEFFICIENT,AS THE FOLLOWING EXAMPLE. 


42 

to (1S) eX 2tie) 
se aa Be Gl Eee [eas (Pa (OPS 2) 
+251 01572) *x%2(072) 
toe ON 2) os) 
Ho" RICO 25" 620072) 
2° XE C37 3) Actes) 
+3™X1 0072) xX 200, 0) 
+U* x) (17 2)" A2(0,2) 
tlt xT C074) i272) 
PUPS 2 a 2) 
tT*eXECL, 207 X2007 1); 


X1(0,3)X2(0,3); 

( input data file; tt5 mvia) 

( Don't care term is 4th term from above with its coefficient '?') 

HOW TO PERFORM LOGIC MINIMIZATION AND SELECT OUTPUT OPTIONS TO GET THE 
DESIRED RESULTS. 

GETTING HELP: THERE IS A HELP ROUTINE FOR THE USERS WHO ARE NOT FAMILIAR 
WITH ALL OPTIONS AVAILABLE IN THIS CAD. 

TYPE THE FOLLOWING COMMAND FOR HELP 


.anviie = ’ 
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(The following HELP messages appears on the screen...... ) 


MVL V1.0 Copyright 1988 Naval Postgraduate School 


usage: mvl [-eimrsvqDPG] [-0.0] [-oFile] [-xFile] [file ...] 
-e - Print the original expression, as parsed 
=f - Print each implicant found 
—m - Print the Karnaugh map of the original and each 


successive expression 
WARNING: impractical for nvar > 3 


Si - Verify the minimal solution 

=§ - Print statistics on selected heuristics 

=a - Quiet, don't print final results 

=0 50 - Build source expressions for all solutions 


whose ratio to the original number of terms 
exceeds this number and output them to a file 
(default is "x.mvl") 


=|) - Execttte the Dueck & Miller heuristic (default) 

=p ~ Execute the Pomper & Armstrong heuristic 

-G —- Execttte the Gold heuristic 

~xFile ~ Output source expressions built using -0.0 to "File" 
instead of "x.mvl" 

-oFile ~ Ouput formatted solutions to "File" 


(Note; In the above HELP messages "print" means "output to CRT".) 


(EXAMPLE 1) USING POMPER AND ARMSTRONG ALGORITHM ON tt2.mvl, OUTPUT ORIGINAL 
EXPRESSION AS PARSED AND FINAL RESULT. 


(command line is as follows...... ) 
ece: /work/mag 
%* mvl -Pe tt2.mvl 


( the following appears on the screen...... ) 


MVL V1.0 Copyright 1988 Naval Postgraduate School 


Expression: 


radix: 4 
nvars: 2 
nterms: 11 
coeff: 2 
AECL 3) 
X2(1,1) 
coeff: 1 
Ki) 
X2(1,3) 
coeff: 2 
x1 2) 
X2(0,2) 
coeff: 1 
Role(s) 
X2(3,3) 
coeff: 2 
X1(0,2) 
X2(0,2) 
coeff: 2 ’ 
X1(3,3) 
X2(2,3) 
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coeff: 3 
X1(0,2) 
X2(0,0) 

coeff: 1 
X1(1,2) 
X2(9, 2) 

coeff: 1 
X1(0,3) 
view 

coeff: 1 
X12) 
X20 lc 

coeff: 1 
X1(1,2) 
X 26071) 


P&A: 37 11 fiplicaitse= 0.77, 
( the last line shows that the POMPER AND ARMSTRONG heuristic produced a 


realization of 3 product terms starting with 11 product terms, fora 
reduction of 27%) 


(EXAMPLE 2) USING DUECK AND MILLER ALGORITHM ON tt2.mvl, OUTPUT 
EACH INPLICANT FOUND AND FINAL RESULT. 
(command line is as follows ... ) 


ece: /work/mag 
So omvl =pi7 ceZ. mull 


#4 
( the following appears on the screen ..... ) 
MVL V1.0 Copyright 1988 Naval Postgraduate School 


D&M MIM: 2*X1( 1, 1)*X2( 3, 3) 
imp: 2*Xi¢ 1, 2)*x2¢ea 3) 


DEM Mite 3*X1C 3, 3)? X20 373) 
imp: 32%xXle 3, 3) x20 ,P -3) 


D&M NIM: 3*X1( 2, 2)#X2( 0, 0) 
Imp: 3*xXLG O.. 2)*% 0 0, 2) 


D&M: Alok implicants - 0.27 


(EXAMPLE 3) USING GOLD ALGORITHM, OUTPUT THE KARNAUGH MAP OF THE ORIGINAL AND 
EACH SUCCESSIVE EXPRESSION AND FINAL RESULT. 
(command line is as follows ...) 
ece: /work/mag 
% mvl -Gm tt2.mvl 
( the following appears on the screen..... =) 
MVL V1.0 Copyright 1988 Naval Postgraduate School 
a ’ 


Orig map (D&M): 
34. 3, “Sa 200 
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a2 3 5 ees | 
S255 3 Sd 
0 Z 2 3 
5. Somes 7 O 
3.3 ee 
33 eo eee 
0 0 0 3% 
3 33. 3 
> a es me: ne 
38 3s ese 4 
0 0 0 4, 
4 eee (0) 
4 AO Ae a 
4 Ae 4s eee a 
0 0 0 4. 
Orig map (P&A): 
Sa 3 EO 
by ee! Ss 
a5 <3 c ae 
0 2 2 Ss. ; 
4. 4 4. O 
4. 4 7 Lake ae 
4. 4 S25 3. 
0 2 2 3 
4. 4 4. 0 
a 4a. Sout om 
4. 4 4. 4. 
0 2 2 4. 
4. 4 4. 0 
4, 4 4A, 4. 
4, 4 , 
0 0 0 4. 
Gold(===): 3/11 implicants - 0.27 


(Note; .(dot) after a number denotes a highest logic value or don't 
care generated in the course of a minimization. Don't cares in the 
original function have no dot.) 


(EXAMPLE 4) USING POMPER AND ARMSTRONG ON tt2.mvl1, OUTPUT THE ORIGINAL 
EXPRESSION AND K-MAP OF THE ORIGINAL EXPERSSION AND EACH 
SUCCESSIVE EXPRESSION WITH FINAL RESULT. 

( commnad line is as follows ...... ) 

ece: /work/mag 

%$ mvl -Pei tt2.mvl 


( the following appears on the screen..... ) 


MVL V1.0 Copyright 1988 Naval Postgraduate School 


Expression: 


radix: 4 
nvars: 2 
nterins: ll 
coeff: 2 
>a res | 
Moet 
coeff: 1 
A173) 
x2C1, 2) 
coeff: 2 
X1(1,2) 
R20, 2) 
coeff: 1 
XE hE. 2) 
X2( 3,3) 
coeff: 2 
ECO) 2 
M2(0 2) 
coeff: 2 
ACS oe) 
M2023) 
coeff: 3 
X10;2) 
X2(0,0) 
coeff: 1 
EE) 
x 2(0e2) 
coeff: 1 
X1(0,3) 
X2015 2) 
coeff: 1 ie 
X1(1,2) 
x01 2) 
coeff: 1 
X1(1,2) 
X2(0,1) 


) 
) 


PG&A MIM: 3*XJ1( 1, 1)*X2( 0, O) 
Imp: 3*X1( 0, 2)*X2( 0, 2) 


P&A MIM: 3*X1( 3, 3)*X2( 3, 3) 
Imp: 34310 3,43)472A2¢ 3) 


P&A MIM: 2*XJ( 1, 1)*X2( 3, 3) 
Imp: 2*X1( 1, 2)*X2( 0, 3) 


P&A: 87 Vi implicants - 0.27 


(EXAMPLE 5) USING DUECK AND MILLER ON tt2.mvl, OUTPUT THE ORIGINAL 
EXPRESSION AND K-MAP WITH FINAL RESULT. 

(command line is as follows...... ) 

ece: /work/mag 

% mvl -Dem tt2.mvl 


(the following appears on the, screen....... ) 


MVL V1.0 Copyright 1988 Naval Postgraduate School 
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Expression: 


radix: 4 
nvars: 2 
nterms: ll 
coeff: 2 
X1(1,3) 
X2(1,1) 
coeff: 1 
X1(1,3) 
RoC 39 
coeff: 2 
ACL 2) 
X2(0,2) 
coeff: 1 
X1(1,2) 
X2(3,3) 
coeff: 2 
X1(0,2) 
X2(0,2) 
coeff: 2 
Al 333) 
X2(2,3) 
coeff: 3 
X1(0,2) 
X2(070) 
coeff: 1 
MEG) 
x210, 2) 
coeff: 1 
X1(0,3) io 
X2(1,2) 
coeff: 1 
X1(1,2) 
X2(1,2) 
coeff: 1 
SLC1G 2) 
X2(0,1) 


Orig map (D6M 
ae, a 


OW Ww 
NG Gy 
I . 


CO Ww Ww Ww 
SO Ww Ww Ww 
Www © 


Oo Ww W& Ww SO Ww Ww Ww 


SO Ww Ww Ww 
Ow Ww Ww 
2» & © 


2 & 2O 


Ona 
Ona. 
Onan 


D&M: 3/11 implicants - 0.27 


8 
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(EXAMPLE 6) USING GOLD ALGORITHM ON tt2.mvl, OUTPUT THE ORIGINAL EXPRESSION, 
K-MAP AND EACH IMPLICANT FOUND WITH FINAL RESULT. 


( command line is as follows...... ) 
ece: /work/mag 
% mvl -Geim tt2.mvl 


( The following appears on the screen....... ) 


MVL V1.0 Copyright 1988 Naval Postgraduate School 


Expression: 


radix: 4 
nvars: 2 
nterms: 11 
coeff: 2 
X1(1,3) 
X2(1,1) 
coeff: 1 
X1(1,3). 
X2(1,3) 
coeff: 2 
RPC) 
X2(0,2) 
coeff: 1 
X1(1,2) 
X2(3,3) 
coeff: 2 
X1(0,2) 
X2(0,2) 
coeffi? 2 
X1(3,3) 
X2(2,3) 
coeff: 3 
X1(0,2) 
X2(0,0) 
coeff: 1 
X1(1,2) 
X2(0,2) 
coeff: 1 
X¥1(0,3) 
X2(1,2) 
coeff: 1 
ALC), 2) 
RZCIG 2) 
eoefft-s1 
X1(1,2) 
X2(0,1) 


Orig map (D&M): 

; ees 8 
a5 
ce 
a 


OW Ww Ww 
NO Ww ty Ww 
Ny Ww Wy 


D&M MIM: 2*X1( 1, 1)*X2( 3, 3) 
Imp: 2+*X1( 1, 2)*X2( 3, 3) 


fad td Ga OE 


Od ad 
Ow le GJ 
OO WwW ly Ww 
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D&M NIM: 34*X1( 3, 3)*x2( 
Imp: 3*X1( 3, 3)#*X2( 


bm 
= 

as 
— 


ors 
3h 
ore 
0 


0 


OW Ww 


4. 
4. 
4 


OW Www 


D&M MIM: 3*X1( 2, Zo x2 
Imp: 3*X1( 0, 2)*x2( 


oOo 


2) 


0 
4. 
4 
4 


© & & 2 
CO & >] & 
CO & & & 


Expression: 


radix: 4 
nvars: 2 
nterms: 11 
coeff: 2 
ACly,.3) 
X2(1,1) 
coeff: 1 
X1(1,3) 
X2(1,3) 
coeff: 2 
AL(1;2) 
X2(0, 2) 
coeff: 1 
eC, 2) 
X2(3; 3) 
coeff: 2 
X1(0,2) 
X2(0,2) 
ecoeris 2 
Ri, 3) 
A702, 3) 
coeff: 3 
X1(0,2) 
X2(0,0) 
coeff: 1 
X1(1,2) 
AZO) 2 ) 
coeff: 1 
X1(0,3) 
AZOL, 2) 
coeff: ] 
X1(1,2) 
wel 2) 
coeff: 1 
ALCL, 2) 
X2(0,1) 
Orig map (P&A): 
Pos oe ome 9 
35 Mics 32 3 
3. wae 3. 2B 
0 2 Ps 3 
P&A MIM: Sex (1, LyX OG, 6) 
Imp: 3*xX1( 0, 2)*X2( 0, 2) 
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eee aa eens 

Ae tA ple, oe 

ey 5 eet ner 

QO “HD 2) 5238 

P&A MIM: 3*X1( 3, 3)*X2( 3, 3) 
Imp: 3*X1( 3, 3)*X2( 1, 3) 

428 A aa 

ee or eee 

4204... “Alea 

O° <27 a 4: 


P&A MIM: 2*X1( 1, 1)*X2¢ 3, 3) 
Imp: 2**1( 1,2) <%2020,s) 


4. 4 4. 0 
4. 4 4. 4. 
4. 4 4 A 
0 0 0 4. 
Gold(===): 2 implicants - 0.27 


(Note; Gold(===) means that both F&A and D&M have the same result and 
the resulting expression of D&M is chosen on ties.) 


(EXAMPLE 7) USING GOLD ON tt2.mvl, OUTPUT STATISTICS ON SELECTED ALGORITHM 
WITH FINAL RESULT. 

(command line is as follows....... ) 

ece: /work/mag 

% mvl -Gs tt2.mvl 

(The following appears on the screen...... ) 

MVL V1.0 Copyright 1988 Naval Postgraduate School 

Gold(===): a7 it implicants ~- 0.27 


Statistics for D&M: 


Num Num Num Avg Num 
Radix Var Expr Terms Term/Expr 
4 2 i 11 11.00 
Eval Comp Pick Comp Gen Next Valid 
Expr Cr MIM RBC Bounds Impl Imp 1 
TOC: 1254 173 4 23 3 Pe! 24 
Avg/Expr: 1254.00 173.00 4.00 Pai ON0) 300 27.00 24.00 
Avg/Term: 114.00 ISt73 0236 2.039 0.27 2.45 2:18 


Statistics for P&A: 


Num Num Num Avg Num 
Radix Var Expr Terms Term/Expr 
4 2 1 11 11.00 
Eval Pick Gen Next Valid 
Expr MIM Bounds Impl Impl 
Tot: a4 4, 3] 30. Py 
Avg/Expr: 534.200 4.00 3.00 30209 27200 
Avg/Term: 48.55 0.36 0-27 273 aaa 
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(EXAMPLE 8) USING GOLD ON tt4.mvl, VERIFY THE MINIMAL SOLUTION BY TEST 
VECTORS SPECIFIED IN THE INPUT FILE. 


(command line is as follows...... ) 
ece: /work/mag 
% mvl -Gv tt4.mvl 


( the following appears on the screen...... ) 

MVL V1.0 Copyright 1988 Naval Postgraduate School 
Gold(===): 3/11 implicants - 0.27 

Verification: 


mC 0,33) 
X2(0,3) 


(0,0) 
(1,0) 
(2,0) 
(3,0) 
(0,1) 
(1,1) 
(2,1) 
(3,1) 
(0,2) 
(1,2) 
(2,2) 
(3,2) 
(0,3) 
(1,3) 
(2,3) 
(3,3) 


uadet a#uadt t @ teu kt &@ 2 


WNN OWWWW WWW WwW OW WwW 


(EXAMPLE 9) USING GOLD ON tt5.mv1l,OUTPUT THE ORIGINAL EXPRESSION, IMPLICANT 


FOUND, K-MAP, STATISTICS ON THE SELECTED ALGORITHM,AND VERIFI- 


CATION OF FINAL RESULT. 

INPUT EXPERSSION HAS A DON'T CARE TERM IN IT. 
(command line is as follows...... ) 
ece: /work/mag 
%$ mvl -Geimsv tt5.mvl 


(The following appears on the screen...... ) 


MVL V1.0 Copyright 1988 Naval Postgraduate School 


Expression: 
a ‘ 
radix: 4 
nvars: 2 
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nterms: li 
coeff: 2 
x, 3) 
X21 1D 
coeff: 1 
XW 
X2(1,;3) 
coeff: 2 
XI(1; 2) 
X2(0,2) 
coeff: 4 
X1(1,2) 
X2(3,3) 
coefr: 2 
ECG; 2) 
X2(0,2) 
coeli=: =2 
ACS 3s) 
X2t2,3) 
coeff: 3 
X1(0, 2) 
X2(0,0) 
coeff: 1 
ALC 2) 
X2(0, 2) 
coeff: 1 
X1CG73) 
X20, 2) 
coeffi. <1 
Adil, 2) 
X22) 
coeff: 1 
AT 2) 
X2(0,1) 





Orig map (D&M): 
5 eee 0 


a Ww W W 


Soe or 
Sa. oh Je 
0 4 a. 


D&M MIM: 3*X1( 3, 3)*xX2( 1, 1) 
Imp: 3*X1( 0, 3)*X2( 1, 2) 


ae 
de De oe 
Fane 
W a & Oo 


D&M MIM: 3*X1( 3, 3)*X2( 3, 3) 
Imp: 3*X1( 3, 3)*X2( 3, 3) 


3 Jeo 0 
4 4. 4 cle 
4 4. 4 ai. 
0 4 4 4 


D&M MIM: 3*X1( 2, 2)*X2( 0, 0) 
Imp: 3+*X1( 0, 2)*x2( 0, 0) 


he de de 


0 
4. 
4 
4 


Oo 2 & 
e e e 
le oe Dp oe 
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Expression: 


radix: 4 
nvars: Z 
nterms: 11 
coeff: 2 
Mitt) 
X2( 1-1) 
coeff: 1 
>A Ge Gat ies 
X¥2(1,3) 
coeff: 2 
AC 12D 
X2(0,2) 
coeff: 4 
X1 (1,2) 
X 23s) 
coeff: 2 
X1(0,2) 
X2(0,2) 
coeff: 2 
K1C37 3) 
X2(2,3) 
coeff: 3 
X1( 0,2) 
X2(0,0) 
coeff: 1 
X11, 2) 
x20C,2) 
coeff: 1 
X1(0,3) 
m2( 1,2) 
coeff: 1 
X1(1,2) 
X2(1,2) 
coeff: 1 
RUC 2) 
X2CO; 1) 


Orig map (P6&A): 
33 von 
3., 73. 
3. Bor 
0 4 


de G) iW 


ye 
0 
2: 
3. 
3. 
P&A MIM: 3*X1( 2, 2)*X2( 1, 1) 

Imp: 3*x1( 0, 2)*xX2( 0, 2) 


we» De te 
Wd GJ ly © 


we De dm 
e e e 


4 
4. 
4. 
0 
PGA MIM: 3*X1( 3, 3)*X2( 2, 2) 

Imp: 3*X1( 1, 3)*X2( 1, 3) 


0 
4. 
4 
4 


Oxbo 
e e ® 
bn Db b & 
Db Db Db & 


Gold(P6&A): 3/11 implicants - 0.27 
Verification: t 


X1(0,3) 
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X2(0,3) 


(0,0) 
(1,0) 
(2,0) 
(3,0) 
(0,1) 
(1,1) 
C271) 
(3,1) 
(0,2) 
(1,2) 
(2,2) 
(372) 
(G53) 
(1,3) 
(2,3) 
(3,3) 


WwW OWWwWwWwW ww Ww Ow ww 


on gd tar aeak uw ht dG 


Statistics for D&M: 


Num Num Num Avg Num 
Radix Var _ Expr Terms Term/Expr 
4 2 1 VBL 11.00 
Eval , Comp Pick Comp Gen Next Valid 
Expr CE MIM RBC Bounds Impl Impl 
Tot: P3344 169 4 31 3 36 33 
Avg/Expr: 1334.00 169.00 4.00 31.00 3.00 36.00 33200 
Avg/Term: $21.27 155.36 0.36 2.82 On27, 3227 3.00 


Statistics for P&A: 


Num Num Num Avg Num 
Radix Var Expr Terms Term/Expr 
4 2 1 11 11.00 
Eval Piek Gen Next Valid 
Expr MIM Bounds Impl Impl 
TOE: git 5 2 94 D2 
Avg/Expr: 811.00 3.00 2200 D4. 00 52200 
Avg/Term: 133 Oey 0.18 4.98 4.73 


(EXAMPLE 10) USING GOLD ON tt2.mvl, IF USER DOESN'T NEED FINAL RESULT 
SPECIFY -q FLAG. 

(command line is as follows....... ) 

ece: /work/mag 

% mvl -Gg tt2.mvl 


( the following appears on the screen...... ) 


MVL V1.0 Copyright 1988 Naval Postgraduate School 


’ 
(ad) HOW TO GENERATE LAYOUT FILE IN MAGIC FORMAT. 
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SPECIFY INPUT , OUTPUT FILE AND THE SELECTED MINIMIZATION ALGORITHM. 
THE FINAL RESULT WILL BE STORED IN THE OUTPUT FILE WITH FORMATTED 


SOLUTION. 


(command line is as follows. input file: tt2.mvl, outfile: tt2.out 
and output file should have -o flag.) 

ece: /work/mag 

%$ mv =G tt2.mvl -obb2zncur 


( The formatted final solution will be saved in an output file which 
will be used for layout generation. If the minimization result is 


worse than the original expression, the original expression will be 
output for producing the layout.) 


( the following appears on the screen...... ) 
MVL V1.0 Copyright 1988 Naval Postgraduate School 


Gold(+==): 37 imolicances — 0.27 


(The formatted output is as follows.) 
# of input f of output 
coeff 
lower upper 
lower upper 
99 (sentinel) 

/* # of input = f of variables 

f# of output = # of expressions */ 


(an example of the formatted output) 


1 


WOO Wh WwW Ws Ww PAD ND 
lad 


9 
( The following command initiates the layout generation 
by using output file ,called tt2.out ) 


ece: /work/mag 
% mvpla tt2.out 


(The following appears on the screen ..... ) 


FRI I II KIRK KKK IKI KHER KKKKEAKEKRKR KKK KKK RK 


Input fetes SOUL 
Output file : tt2.out.mag , 


RHI IKI RIKI KKK KKK KKKKKKKKRKKKRKKR KKK KKK KKK 
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Check: 
2 inputs 1 outputs function 


fOr 2x0(1 72) x1C373) 
2 8x03 73) x11; 3)) 
+ 3x0(0;,2)x1(G0, 2) 


Now, generating CMOS-MV-PLA layout in magic format...... 


£0", 2X0 (1, 2) X11 378) 
PreaxO(3 3) x1 (173) 
* 3X0 (07 2)x1(052) 
Output to Gt2 out. magus. 


Done ! 


hhakhhhkhhteh EXIT ttt kt tthe kik 


( Once layout file with the Magic format is created, 
the user needs to invoke the Magic graphic package. ) 

( command line is as follows..... ) 

ece: /work/mag 

% Magic tt2.out 


(The following appears on the screen eee ) 
Magic - Version 4.10 - Last updated 12/6/85 at 18:33:06 
Using technology "scmos". 


£t2.out: 500 rects 


(After the whole layout is displayed on the CRT, the user can create a file 
for use in generating a hard copy of the layout and exit by typing the 
the following command to create the .cif file for plotting.) 


:c1if write tt2.out 
>:quit 


(The following command line allows the user to have a hard copy of layout.) 
ece: /work/mag 

SeeCiiplot —r LtZ outer 

{1} 20267 


ece: /work/mag 

% Window: -10950 12600 -2100 51000 

Scale: 1 micron is 0.0198588 inches (504x) 
The plot will be 0.39 feet 


{1] Done cifplot -I -P /tools/berk86/lib/patterns -r tt2.out. 


(the following is an example of a hard copy of MVL-PLA layout. ) 
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 cnenaedetedeendods- anadkedecadh arakemdcantvandeaamds dead deaetaaiaden dead-end 


(RUTRIEUUPLEUOUICS USE IE Pee eee ees eee eee ee eee 
Ses RAE Pine 1 PIF ARSE H eS Sm SS ine te OME ' 
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_ "ett 
TAODUIA 23015415 


BEIZ1 SS6Z1- 


BBlc- 


BBSB6IB°S $) YOuDjm | 20, e>$ --- gBBls 


wma KK KK KK KK KK KKK KKK KKK KKK KKK AAR KRHA HHAHAKKHaEKHKKaEKaR KKH Ha KKK KKK KKK 
Bem KK KH KK KKK KI KKK IKK KKK KKK KKK KKK KKK KKK KKK HI K KKK KKK KKK KKKKKEKKEKEKKEARKEKKARS 
FRR RK RK KR KK KK KK RIK RK RK KKK KK KHIR KKK KK HK KKK RK KIRK KKK RERKKKRKKKKEKEKKRKAKKRKRKR 


THIS COMPLETES THE DESIGN PROCEDURE FOR THIS CAD TOOL. 


RK KK IK RR RK RK RK KK RR KKK KK KK KK KK KKK KKK RK KKK HR RIK RRR KKK KKEKHAKAKAKAKHKKK a 
KR RR HK KH KR KK KKK KKK KKK KKK KK HK KK KKK RK RK KKK RK KHER RK RK RKKAKRKERKAKRKAKKKKARHKKaKRKaKRaK Ka 
KKK RR RK KKK KKK KIRK RK KR KKIKKRAKRKRKRKKKKKEKEKEKKKKKKKRKRKKKAKRaAKKKKKKKKRE 


2.ANALYSIS PROCEDURES 
THERE IS A RANDOM EXPRESSION GENERATOR IN THIS CAD. TIE USER CAN GENERATE 


A SET OF RANDOM EXPRESSIONS OF A SPECIFIC CLASS AND ANALYZE THE DIFFERENT 


MINIMIZATION PROGRAMS. 


40) 





(A) HOW TO GENERATE RANDOM EXPRESSIONS 


(EXAMPLE 1) THERE IS HiELP MANUAL FOR USERS WHICH CAN BE OBTAINED AS FOLLOWS. 


ece: /work/mag 
% mvltest -- 


(The following HiELP messages appears on the screen..... ) 
MVLTEST V1.0 Copyright 1988 Naval Postgraduate School 


usage: mvltest [-sN -rN -iN -tN -oN] 


-sN - Cycle the random seed N times (default = 0) 
—rN - Set radix to N (default = 4) 

=aN - Set number of inputs to N (default = 2) 

= tN - Set number of terms to N (default = 1) 

—ON - Set number of output to N (default = 1) 


(EXAMPLE 2) THE FOLLOWING COMMAND INITIATES A RANDOM FUNCTION GENERATOR 
TO CREATE A SET OF RANDOM EXPRESSIONS. IT CYCLES THE RANDOM SEED 50 
TIMES AND GENERATE 20 EXPRESSIONS OF 4-VALUED 2-VARIABLE LOGIC EXP- 
RESSION. ALL OF THEM ARE STORED IN tt6.mvl AND HAVE 10 PRODUCT TERMS. 
THE NUMBER OF OUTPUTS IS THE NUMBER OF FUNCTIONS IN THE SAMPLE SET. 


(command line is as follows.!'.. ) 
ece: /work/mag 
% mvltest -s50 -r4 -i2 -t10 -020 > tt6.mvl 


(The following is a set of random expressions generated by 
the previous comnand. ) 


4:2: 
+2*X1(1,3)*X2(3, 3) 
Hd x1(0,0)*X2(1, 2) 
oie X0(0,3)*X2(1,2) 
Freee C1, 2) *X2(1, 2) 
toe AI (1,2) *X2(0,)) 
+2*X1(0,0)*X2(1,1) 
Toe xO, 1) *X2(1,3) 
+2*X1(3,3)*X2(0,1) 
+3*X1(0,3)*X2(1,1) 
+2*X1(0,3)*X2(0,0); 


+41*X1(0,3)*X2(0,3) 
#2*X1(1,1)*X%2(0,3) 
+1*X1(1,2)*X2(0,1) 
eU*X1(1,1)*X2(0,2) 
feo x1 (0,0)*X2(2, 2) 
fee X15 3,5) *X2(0,3) 
4#2*X1(1,2)*X2(1,1) 
#2*X101,1)*X2(2,2) 
+9*X1(3,3)*X2(0,3) 
#U*X1(1,1)*X2(1,1); 
4:2: 

*2*X1(1,2)*X%2(0,1) "3 ; 
mo PCT, 2) *X2(0, 3) 
pot Al(2,3)*X2(0,1) 


4] 


+2*XT CO, 1) XD oe 
+2*X1(1,3)*X2(2,2) 
#14X1(0,2)*X2(3, 3) 
+3*X1( 1, dpe xo orn 
+1*X1(1,2)*X2(0,0) 
+2*X1(0,0)*X2(0,3) 
+3*X1(3,3)*X2(0,2); 


+2*X1(1,2)*X2(3,3) 
+2*X1(1,2)*X2(2,3) 
+1 4X1 00, 14 X201 on 
+1*X1(1,3)*X2(3,3) 
+3*X1(0,3)*X2(1,1) 
#1*X1(3,3)*X2(2,3) 
+3*X1(3,3)*X2(1,2) 
+P* S17 (1,0) 2% 20570) 
+2*#X1(0,1)*X2(3, 3) 
#3°X100,1)* x 201306 


42° XV 3), 3)? Xeon) 
+3*X1(2,2)*X2(3,3) 
#1*X100, 1)* X20 
+9*°X101, 3) *X202 oy 
+2*X1(1,1)*X2(0, 2) 
+3*X1(0,2) *X2(2,3) 
+1*X1(2,2)*X2(1,1) 
+T*X1 (2,2) hol, ) 
$1*XE(2, 2) oot 
+3*X1(0,1)*X2(1,2); 


#3* X10 152) 4201 2) 
43° X1 (2,9 )*X2¢ioy 
+1*X1(2,2)*X2(0,0) 
+3*X1(1,2)*X2(1,1) 
#1*X1(3,3)*X2(1,2) 
+3*X1(1,3)*X2(0,0) 
+2*X1(0,2)*#X2(1,2) 
+e XT COnd x20 ee) 
+1*X1(0,1)*X2(0,0) 
+A X13) 3X21 0710) 


+3*X1(0,3)*X2(1,3) 
+3*X1(1,1)*X2(1,3) 
+2*X1(0,3)*X2(0,2) 
+3*X1(2,3)*X2(0,1) 
+2*X1(0,1)*X2(3,3) 
+3*X1(0,2)*X2(3,3) 
+9* X11), 1) *X 200, 3) 
#3*X1(1,1)? X20, 1) 
429X110 1) 2X2 iy 
+deXI12, 2)? XO as 


+3*X1¢1,2)*X2¢1,2) 
4+14*X1(0,3)*X2(2,3) 
+24#X1(3,3)*X2(2,3) 
+3*X1(0;3) 742015) 
+2*X1(2,3)*X2(2,3) 
+2*X1(0,0)*X2(0,0) 
+#3*X1(0,0)*x%2(1,1) 
+3*X1(0,1)*X2(1,2) 
+34*X1(2,3)*X2(0,0) 
+3*X1(2,3)*X2(0,1),; 
4:2: 
+34*X1(2,3)*X2(0,3) 
+3*XU0C 2) X27 (052) 
+3*X1(1,2)*X2(2,2) 
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#1 *X1(0,0)*X2(2,3) 
+#1*X1(0,3)*X2(2,2) 
4$2*xX1(1,3)*X2(0,0) 
FOS Xe eX 2 1, 3) 
+14*X1(0,2)*X2(0,1) 
+3*X1(0,1)*X%2(3,3) 
£2 epee Gl) 7 
4:2: 
+3*X1(2,2)*X2(0,1) 
+14*X1(0,2) *X2( 2,2) 
#3*X1(0,3)*X2(2,3) 
#2*X1(0,1)*X2 (1,3) 
+2*X1(0,0)*X2(1,2) 
41*X1(3,3)*X2(0, 3) 
2X 3) X2( 2,3) 
+1*X1(0,3)*X2(0,0) 
+2 ©, 3) *X2( 1,3) 
+35 *X CC, O)*xX2(1,1); 


+1? MgO 3 * X2 (2,2) 
+3*X1(1,3)*X2(0,3) 
+1* 8102, 2)*X2(0, 2) 
+1*X1(1,1)*X2(0,1) 
41*X1(1,3)*X2(1,1) 
+1*X1(3,3)*X2(2,3) 
+3*X1(2,3)*X2(2,3) 
E2eX NCO) *xX2¢1,1) 
+14*X1(1,3)*X2(0, 2) 
+14*X1(0,0)*X2(3,3); 


tC ye X20270)) 
+ 2% M3) *X2(0,2) 
+3* X11, 3) *X2(2,3) 
+1*X1(1,3)*X2(2,3) 
4#1*X1(0,2)*X2(0,1) 
+3* XV 007) *xX2(0;,0) 
EIN M ee eel 
+3* 300252) *X2(0, 3) 
le xv 3) 7A2(0,0) 
FV * SACO A *X 200, 2) 4 
4:2: 
Ue X UO 3X 20359) 
+2*XU( 073) * X20 272) 
+1 * X07 3) *X2(073) 
USAW) x2 (07) 2) 
+1* X23) Xx 2(0,3) 
+I* XPT 3) *X2(0, 2) 
+2* X00, 1) X20, 2) 
to X22) X23, 39) 
#1*X1(1,1)*xX2(0,3) 
+3*X 101, 207X202, 2)3 


+2*X1(0,0)*X2(0,0) 
PLY XT 273) x 20173) 
+3*X1(2,3)*X2(0,3) 
+1*X 133) X 20073) 
+2*XU(373)*k2(17 2) 
+3* X10, CO) 4X2 (252) 
+2*X1(0,1)*X2(0, 2) 
+2* X02 3) Xe 173) 
+2*X UG X21, 2) 
+2*X1(0,1)*X2(3,3); 
4.22: 
Perl 3)t x20) 71) 
+2*X1(0,2)*X2(0,3) 
+#1*X1(0,3) *X2(3,3) 
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£2*X1(0,3) *x2(070) 
+3*X%1(0,2)*X2(0,3) 
+1*#X1(1,3)*X2(1,2) 
+1*X1(0,0)*X2(1,3) 
+3*X1(2,2)*X2(0,0) 
43*X1 02,3) * x ota 2 
+2*X1(0,3)*X2(0,2); 


+2*XUC 253) * RoC ee) 
+2*R1 (1,0) * 5 ot 
+3*x1¢1, 3)? Xoct 
+24#X1(3,3)*X2¢1,1) 
+1*X101, 1) *X203, 2) 
43*XI(0,3)*X2(0,1) 
+34N1 02, 3) R213) 
41*X9(39/ 3)? X21) 
re (lp ys) rele 25 
+IERM 22) PRON Dae 


+3*X1(1,2)*X2(0,0) 
+34 5070, 1) *#X2(37 9) 
+3912, 3) * X21 22> 
42% MIN 33) 4x23 795 
+1*X1 (009% X2(1,3) 
+14 X10) 2 22) 
414X100; 3) #201) 
4+34*X1(0,1)* 20D) 
+3*X1 00,3) X201 2) 
#2*X100,3)4*Xo(27 


43*X1(0,3)*X2(0,0) 
Pye xo) 2 
+3*X1(0,1)*X2(2,2) 
+1*X102,2)*%2(1, 1) 
+2*X1(373 "2 (ed 
+3*X1(0,3)*X2(1,2) 
$1 *X1 Cie Xe ee 
tO PRU Ue) x 2007) 
+3*X1(1,3)*X2(0,3) 
+14*X1(3,3)*X2(2,3); 


+3*X1(0,3)*X2(0,0) 
+IF NTE, Stet O, 2} 
+3* X10) 3)" x20 ain 
+2*X1(0,1)*X2(1,1) 
+1*X1(0,2)*X2(0,0) 
+1*X101,1)*X2(0,1) 
#2*X1(3,3)*X2(3,3) 
+2#X1(0,3)*X2(3,3) 
+3*X1¢1,3)*X2¢1,1) 
+24X1(0, 2)*X2¢0,0)) 


#3*X1(2,2)*X2(1,3) 
419X170 2,99 9% 2 (oe 
41*X1 (0,0) *X2002) 
+2*X1(07 1) *X2005 3) 
+3*X1(0, 3) 7 X2(07 0) 
+2*XUCL Ves 
+3*X1(0,3)*X2(0, 2) 
+3*X1(2,3)*X2(0,0) 
+1*X1(0,2)*x2(0, 2) 

+2*X1 (1,398 82 00m 
(data file; tt6.mvl) 
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U 
(EXAMPLE 3) THE FOLLOWING COMMAND LINE APPLIES Gold ON tt6.mvl. THE 
STATISTICS AND THE FINAL RESULTS ARE FOLLOWED. 


ece: /work/mag 
% mvl -Gs tt6.mvl 


MVL V1.0 Copyright 1988 Naval Postgraduate School 


Gold(P&A): 6/10 implicants - 0.60 
Gold(===): 6/10 implicants - 0.60 
Gold(===): 3/10 implicants - 0.30 
Gold(===): 37 1G implicants - 0.50 
Gold(===): 4/10 implicants -— 0.40 
Gold(D&M): 4/10 implicants - 0.40 
Gold(===): 3/10 implicants - 0.30 
Gold(===): 4/10 implicants - 0.40 
Gold(===))+ 1/10 implicants =O. 10 
Gold( ===): 7/10 implicants — 0.70 
Gold(===): 3/10 implicants - 0.30 
Gold(PS&A): 6/10 implicants - 0.60 
Gold(===): 3/10 implicants - 0.30 
Gold(===): 5/10 implicants - 0.50 
Gold(===): 27 10 implicants —- 0.30 
Gold(===): 4/10 implicants - 0.40 
Gold(D&M): 4/10 implicants —- 0.40 
Gold(===): 2/19 implicants —- 0.20 


Gold(===): 5/10 implicants —- 0.50 
Gold(===): 3/10 implicants - 0.30 


Statistics for D&M: 


Num Num Num Avg Num 
Radix Var Expr Terms Term/Expr 
4 2 20 200 £0500 
Eval © Comp Pick Comp Gen Next Valid 
Expr CE MIM RBC Bounds Impl Impl 
HOt: 31979 3012 101 943 81 904 970 
Avg/Expr: 1398.95 150.60 5.05 47.15 4.05 45.20 48.50 
Avg/Term: 159.89 13..06 eet 4.71 0.40 4.52 4.85 


Statistics for P&A: 


Num Num Num Avg Num 
Radix Var Expr Terms Term/Expr 
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4 2 20 200 10.00 


Eval PIek Gen Next Valid 
Expr MIM Bounds Impl Impl 
Jot: 16638 101 81 1032 JoO33 
Avg/EXxpr: 831.90 5.05 4.05 51.60 51.65 
Avg/Term: ele ke) Ov5n 0.40 5.16 Sly, 


( EXAMPLE 4) HERE GENERATE ANOTHER SET OF RANDOM EXPRESSIONS, 10 EXPRESSIONS 
OF 4-VALUED 3-VARIABLE LOGIC EXPRESSION AND SAVE THEM IN THE DATA FILE 
tb] -mvi. 

(command line is as follows...... ) 


ece: /work/mag 
% mvitest -r4 -13 -t5 -ol0 > EEF nt 


(EXAMPLE 5) APELY Gold ON tt7.mvl, GET THE FINAL RESULTS AND STATISTICS. 
(commmand line js as follows..... ) 

ece: /work/mag 

% mvl -sG tt7.mvl 


(the following appears on the screen 5 ) 


MVL V1.0 Copyright 1988 Naval Postgraduate School 


Gold(===): 4/7 implicants — O.57 
Gold(D&M): 6/7 impricants = 0.86 
Gold(===): Wad) implicants — 1.00 
Gold(P&A): 8/7 implicants ~ 1.14 
Gold(D&M): ee implicants — 1.00 
Gold(D&M): 5/7 implicants ~ 0.71 
Gold(D&M): 6/7 implicants - 0.86 
Gold(D&M): Td implicants — 1,00 
Gold(===): aft implicants — 1.00 
Gold(===): iva impli eatits ae 100 


Statistics for D&M: 


Num Num Num Avg Num 
Radix Var Expr Terms Term/Expr 
4 3 10 70 7,00 
Eval Comp Pick Comp Gen Next Valid 
Expr Ch MIM RBC Bounds linpL Impl 
Toe 54192 2525 74 1126, 64 922 1257 
Avg/EXpr: 5419.20 Zoe. 0 7.40 Ti2<60 6.40 92.20 125-70 
Avg/Term: gfe ee ee 36.07 106 16.03 Onot hae re Be 17.96 
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Statistics for P&A: 


Num Num Num Avg Nun 
Radix Var Expr Terms Term/Expr 
4 2S 10 70 T2200 
Eval Pick Gen Next 
Expr MIM Bounds Impl 
Tot : 25434 19 69 1269 
Avg/Expr: 2543.40 7.90 6.90 126.90 
Avg/Term: 363.34 ere 0.99 15.15 


Valid 

Impl 
P3333 
33250 

19.04 


(EXAMPLE 6) OUTPUT SOURCE EXPRESSIONS FOR ALL SOLUTIONS WHOSE RATIO TO THE 


ORIGINAL NUMBER OF TERMS EXCEEDS 1.00. ALL WILL BE STORED IN tt7.out 


WHICH HAS -x FLAG ON IT. 
(command line is as follows...... ) 
ece: /work/mag 
% mvl -1.0 —-xtt7.out -Gs tt7.mvl 


(the following appears on the screen...... ) 


MVL V1.0 Copyright 1988 Naval Postgraduate School 


Gold(===): 4/7 implicants — (MS, 

Gold(D&M): 6/7 implicants - 0.86 

Gold(===): 7/7 implicants - 1.00 

Gold(P&A): 8/7 implicants - 1.14 

Gold(D&M): Pr | implicants - 1.00 

Gold(D&M): 577 implicants - 0.71 

Gold(D&M): 6/7 implicants - 0.86 

Gold(D&M): 7/7 implicants - 1.00 

Gold(===): oe implicants - 1.00 

Gold(===): 1m a; implicants - 1.00 
Statistics for D&M: 

Num Num Num Avg Num 
Radix Var Expr Terms Term/Expr 
4 3 10 70 7.00 
Eval Comp Pick Comp 
Expr CE Mil RBC 
Tot: 54192 2520 74 1126 
Avg/Expr: 5419.20 292.90 7.40 112.60 
Avg/Term: 774.17 36.07 1.06 16.09 
Statistics for P&A: = 
Num Num Num Avg Num 


Gen 
Bounds 
64 
6.40 
0.91 


Next 
Impl 
922 
o2. 20 
lt je BY 


Valid 

Iinpl 
257 
125270 

17.96 


Radix Var Expr Terms 
4 3 10 70 
Eval Pick 
Expr MIM 
Tot: 25434 7 
Avg/Expr: 2543.40 7230 
Avg/Term: 363.34 13 


Term/Expr 


Next Valid 
Impl Impl 
1269 1393 


126290 133,50 
U3 a Bs, 19.04 


(EXAMPLE 7) THE FOLLOWING SHOWS ALL THE EXPRESSIONS WHICH EXCEEDS THE 


h FS 
4:3: 


SPECIFIED RATIO IN THE PREVIOUS COMMAND LINE. 


A:8/7 


+2*x1(1,2)*x2(1,2)*x3(0,1) 


F2*x1 (0,2) x21 3 es ee 
+1**1(0 094201) xo 
+3tX1 (1, 3) X23 I 
ee xt (1, 1) tee (0,3) xo 
+TAX 0073) *x 200 73) * x3 Oe) 
t1*x1(1,1)*x2(0,3)* x32(0, 1); 


f D&M: 8/7 


4:3: 


t1*xl 01,1) *x200, 1) *x32 3) 


+2*xXE( 1,3) *x2 2, 2)% x3( 073) 
Wee eee as) a 
) 


+2*x1(0,2)*x2(0,0)*x3(0, 


+3*x1(3,3)*x2(0,2)*x3(0,1) 
+3*x1(3,3)*#x2(1,1)*x3(2,2) 
+2*x1(2,3)*x2(1,1)*x3(1,3); 


fF P&A: 8/7 


4:3: 


+2*x1(0,3)*x2(0,3)*x3(1,1) 


#1*x1(3,3)*x2(2,3)*x3(1,3) 
+3*x1(1,2)*x2(0,1)*x3(2,2) 
+3*x1(1,1)*x2(3,3)*x3(2,2) 
+3*x1(2,2)*x2(3,3)*x3(0,3) 
+2*x1(0,1)*x2(1,3)*x3(1,1); 


A P&A: 8/7 


4:3: 


+3*X1(00, 2) *X2(375)" X30) 
+1*xi (1,1) °%x20073)7 x30, 2) 
+2*x1(1,1)*x2(1,1)*x3(1,2) 
+14x1(0,1)*x2(1,1) 4337) 
#2*x1(0,1)*x2(2,3)7 xo ae 
#2*xX1(0,1)*x2(0,1)*7x300,0) 
+3*X1( 2759 X21, 3) Xo 2) 


(data file; tt7.out) 
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APPENDIX B 


Program Listing 


THE FOLLOWING ARE THE MODULES USED FOR THE MVL-CAD. THE ENTIRE PROGRAM 


IS ABOUT 4,000 LINES LONG. 


HERE IS A BRIEF OVERVIEW OF WHAT IS INCLUDED. 


(1) 


(2) 


(3) 


i) 


(5) 


VERSION BANNER, HWELP DISPLAY AND GLOBAL DATA STRUCTURES FOR MVL-CAD. 


main.c 

THIS MODULE 1) PROCESSES COMMAND LINE ARGUMENTS 2) OPENS FILES AND 
3) CONTROLS THE PARSING, LOGIC MINIMIZATION AND OUTPUT FUNCTIONS. 
ALSO, IT INCLUDES Gold(), WHICH IMPLEMENTS THE GOLD HEURISTIC ( 
WHICH CHOOSES BETWEEN POMPER-ARMSTRONG AND DUECK-MILLER, PREFERRING 


DUECK-MILLER ON TIES. ) 


defs.c 


THIS MODULE INCLUDES ALL GLOBAL VARIABLE DEFINITIONS. 


parse.1] 

THIS MODULE IS A LEX SPECIFICATION FILE FOR THE INPUT SCANNER/PARSER 

FOR THE MVL-CAD. IT SPECIFIES THE GRAMMAR OF EXPRESSIONS INPUT BY THE 
USER. ALSO, IT INCLUDES verify(), WHICH ALLOWS THE USER TO DETERMINE 


THE FUNCTION VALUE STORED. 


dm.c 
THIS MODULE IMPLEMENTS THE DUECK-MILLER HEURISTIC FOR MINIMIZING AN 
EXPRESSION. 
THE FOLLOWING SUBROUTINS ARE INCLUDED IN THIS MODULE. 
a) Dueck_Miller() 
- performs the Dueck and Miller heuristic on the input expression 
and returns the number of implicants used to cover it. 
b) c£(E,X) 
- computes the clustering factors at X, where X is a vector of 
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coordinates of a non-zero minterm which has the smallest value. 
- returns the clustering factor, an integer. 
c) *mim(E) 

- finds the most isolated minterm determined by the minterm with 
smallest CF and returns a vector of integers representing the 
coordinate of it or NULL if no more minterms 

- the function value at that location is also returned as the last 
integer in the vector. 

ad) valid_implicant(1I) 

- checks implicant validity and returns 1 if the proposed impli- 
cant is valid; else returns 0.(An implicant is valid if it is 
indeed an implicant of the function. ) ) 

e) Soihpueeetberaey) 
- accepts implicant I and calculates the relative break count of it. 
£) “pick implicant (x) 

- finds the most isolated minterm according to the smallest CF. 

—- generates all implicants that covers this minterm. 

~- computes the rbc for each and chooses the implicant with the 
smallest rbc (best implicant). 


- returns the best implicant. 


(6) pa.c 
THIS MODULE IMPLEMENTS THE POMPER-ARMSTRONG HEURISTIC FOR MINIMIZING 
AN EXPRESSION. 
THE FOLLOWING SUBROUTINES ARE INCLUDED IN THIS MODULE. 
a) Pomper_Armstrong( ) 
- performs the Pomper-Armstrong heuristic on the input expression 
and returns the number of implicants used to cover it. 
b) *mim(E) 
- finds the most isolated minterm in the input expression based 
on the random weight generated by random number generator. 
~ returns a vector of integers representing the coordinates of 


the most isolated minterm or NULL if no more minterms. 


- the value of the function corresponding to the most isolated 
minterm is also returned as the last integer in the vector. 
Cc) random() 
—- returns a random number using a linear congruential method. 
qd) valid _implicant(1I) 
- checks implicant validity and if valid, returns 1 and the num_ 
zeroed, the number of minterms driven to zero by this implicant 
when subtracted, else returns 0. 
e) *pick_implicant(X) 
- finds the most isolated minterm randomly and then produces all 
implicants that cover this minterm. 
— computes num_covered, the number of minterms covered by this 
implicant. 
- picks the best implicant which drives the maximum number of 
minterms to zero. If there are more than one such implicant 
the largest implicant is chosen. 


- returns the best implicant. 


(7) common.c 


THIS MODULE INCLUDES HEURISTIC SUPPORT FUNCTIONS COMMON TO DUECK-MILLER 
AND POMPER-ARMSTRONG. 
THE FOLLOWING SUBROUTINES ARE INCLUDED IN THIS MODULE. 
a) * eval(EB,X) 
- evaluates the expression at (x1,x2,x3,....,xn) 
b) *next_coord(T) 
- computes the next possible coordinate for term. 
- returns an integer vector(X) containing the coordinates. 
c) *gen_bounds(X) 
- generates the permissible bounds around location X in the work- 
ing expression for use in picking the next implicant that covers 
a 
—- returns a bounds array. 


d) *next_implicant(B) 


ol 


- chooses X as the first implicant. 
- On each call, returns the next implicant with a bound array(B) 
from *gen_bounds() that specifies the implicant. 
e) subtract implicant(1) 
—- adds the best implicant(I) to the working expression as a neg- 
ative term (negated coefficient). 


- Also adds the best implicant to the final expression. 


(8) alloc.c 
THIS MODULE INCLUDES ALL MEMORY ALLOCATION FUNCTIONS, SUGHAAS alice 
term(), alloc_bound(),an1 dealloc(),etc. IT IS USED TO ALLOCATE SPACE 
FOR A TERM ARRAY, BOUND ENTRIES AND DEALLOCATE OR DUPULICATE THE EXP- 


RESSION. 


(9) mvitest.c 
THIS MODULE IS A TEST SOURCE FILE GENERATOR. IT GENERATES ANY TYPE OF 


MVL EXPRESSIONS ACCORDING TO THE SPECIFICATION GIVEN BY THE USER. 


(10) Makefile 
THIS MODULE KEEPS TRACK OF WHICH SOURCE FILES ARE DEPENDENT ON OTHERS 


AND PROVIDES FOR AN EFFICIENT RECOHPILATION OF A PROGRAM. 


jaa] 
bo 


/* mvl 
- A multi-valued logic expression compiler/analyzer/optimizer 


LCDR John M. Yurchak 
LT Hoon Seop Lee 
Dr. Jon T. Butler 


Copyright 1986 U.S. Naval Postgraduate School, Monterey, California 
All rights reserved a / 


finclude "“defs.h" 
finclude <ctype.h> 


/* Version banner and help display */ 
static char 


*version = "“\nHVL V1.0 Copyright 1988 Naval Fostgraduate School\n\n", 
usage = "usage: mvl [-eimrsvqDPG] [-0.0] [-oFile] [-xFile] [file ...J)\n", 


*help{) = { 
-e - Frint the original expression, as parsed", 
. -i ~ Frint each implicant found", 
~m - Frint the Karnaugh map of the original and each", 
Hs successive expression", 
° WARHING: impractical for nvar > 3", 
ig -v - Verify the minimal solution”, 
“ ~3 - Frint statistics on selected heuristics", 
= -qg - Quiet, don’t print final results", 
: = - Permit negative coefficients in terms", 
= -0.0 - Build source expressions for all solutions", 
" whose ratio to the original number of terms", 
“4 exceeds this number and output them to a file", 
" (default is \"x.mvl\")", 
Ww -D - Execute the Dueck & Miller heuristic (default)", 
as -FP - Execute the Fomper & Armstrong heuristic", 
-G - Execute the Gold heuristic", 
5 ~xFile - Output source expressions built using -0.0 to \"File\"", 
" instead of \"x.mv1\"", 
is -oFile - Ouput formatted solutions to \"File\"", 
NULL 


o3 


/* Global data structures --9- eee en rrr rrr rene nn-= ie / 
/* Logic expressions: 


E orig 
- holds the original input expression as parsed 


E work 
- a copy a E orig 
- implicants are subtracted from this expression as terms 
during the coures of optimization 


E final() 
- the result expression (starts out empty) 
- each term is one implicant found during optimization 
- each heuristic has its own E_final (for comparison) 


etd 


Expression 
F orig = [ 0,0,0,NULL }, 
E work = { 0,0,0,NULL }, 
E final[(2) = [ 
{ 0,0, 0, NULL. }, 
( 0,0,0,NULL } 
; 


int HEUR; /* Current heuristic 
* HEUR indexes into E_ final[() 
* depending upon the currently 
* active heuristic 


a 
int FINAL; /* Index of the selected final 
* expression 
a7 
char 
expression(MAX BUF+1), /* Parser working buffer */ 
if name (MAX PATHH1), /* Input file name */ 
xf name(MAX_FATH+1) = "x.mv1\0", /* Expression output file name */ 
*tmp name = “mv lLXXXXxXX", 
of name (MAX_PATHH1); /* Output file name (for result) */ 
int token type; /* Scanner token code */ 
int completed; /* Flag: 
* 1 if a complete sentence has been 
* recognized, 0 otherwise 
ae / 
jmp_buf eof_context; /* Saved context for lex end of file */ 
/* user options */ 
int © flag = 0, 
e tlag = 0, 
i flag = 0, 
m flag = 0, 
s flag = 0, 
v_ flag = 0, 
D flag = 1, 
P flag = 0, 
G flag = 0, 
q flag = 0; 
int nvar, /* Number of variables in the current expression */ 
radix; /* Radix of the current expression */ 


qn 
bbnn 


int tot expr = 0, /* Total number of expressions optimized */ 
tot_terms = 0; /* Total number of terms in all input expressions */ 


MVL_ stats 
*STAT, /* Fointer to the current statistics record */ 
pti stat, /* Stats record for the D&M heuristic */ 
FA_stat; /* Stats record for the PGA heuristic */ 


float FO ratio = 0.0; 
FILE *xfiles: 


/* 
- yywrap() is called by lex on EOF. 
- Thus, we execute a non-local goto back to the setjmp call which 
saved the context, at which point we open the next input file. 
- Notice that yywrap(} never actually returns. 

yh 

yywrap () 

{ 
longjmp (eof _context,1); 


4 


a15) 


main(argc, argv) 
char *argvt); 
[® ~--~--3- === = 8 a a ee 
t: function: 
- Process command line arguments 
- Open files 
- Control the parsing, optimization and output functions 


ee we we ee ee ee ee ee ee ew ee we ew ee ee ee eee ee ee ee ee ee ew we ee ee we ee we 


register ap; 
char *p?; 


fprintf(stderr, "%s",version); 
/* Do the options */ 


for {ap-l; ap < argc; apt+) { 
Pp © argv{ap}; 
Lf [(* ptt a=" —*) 
/* loop through an option string */ 
while (*p) { 
switch (*pt+) { 
case 'D’': 
D flagtt+; 
break; 
cage "F's: 
e f£lagtt; 
4£ (D flag==1) 
Db flaq=-; 
break; 
case 'G': 
G_ flagtt; 
D flagtt+; 
Pflagtt; 
break; 
case ‘vy’: 
vo tlagts? 
break; 
case ‘x’: 
if {*p) { 
strcpy (xf_name,p); 
/* terminate the while loop */ 
Spo eSe Noe 


} 

else { 
fprintf(stderr, "mvl: filename expected after -x\n"); 
fprintf(stderr, "%s",usage); 
exit (2)? 

) 

break; 

case 'o’: 

a =) A 
strepy {of_name,p); 
/* terminate the while loop */ 
ahh, - oe 6 hake 

) 

else { 
fprintf(stderr, "mvl: filename expected after -o\n"); 
fprintf(stderr, "%9", usage) ; 
exit (1); 

} 

break; 

case ‘q's: 
Gg tage, 
break; 


case ‘3s’; 

Ss flagtt; 

break; 
case ‘r’: 

ry flagtt; 

break; 
case ‘m’': 

m_ flagtt; 

break; 
case ‘e’: 

e flagtt; 

break: 
case ‘i's 

i_flagtt; 

break; 
case ’,’ 
case ‘'Q’ 
case ‘1’ 
case ’2' 
case ‘3’ 
case ‘'4’ 
case ’5’ 
case ’ 6’ 
case ‘7° 
case ' 8’ 
case '’9’ 

p--? 

sscanf (p, "Sf", 6FO ratio); 

‘p= "\0'; 

break; 
case ’~'3 

fprintf(stderr, "%s", usage); 

for (ap=0; helplap); apt+) 

fprint£ (stderr, "¢s\n",help{ap)); 

exit (0)? 
default: 

fprintf(stderr, "mvl: unknown switch - %c\n",p[~-1])); 

fprintf(stderr, "$s", usage); 

exit (1); 


ee se eo ee 


ee 


/* Erase the old expression file (if there is one) */ 
unlink (xf_name) ; 


/* If the user specified an output file ... */ 
1£ (of name{0J) |{ 
/* erase the existing output file and open the new one (append) */ 
unlink (of name) ; 
mktemp (tmp name) ; 
if ((yyout=fopen(tmp_name,"a")) == NULL) { 
fprintf(stderr, "mvl: can’t open temp file - %s\n",tmp_name) ; 
exit (l}; 


/* Init the scanner */ 
token_type = TOK_EOL; 


/* For each input file */ 

for (ap-l; ap < argc; aptt) { 
p = argvlap};? 
/* skip the option arguments */ 
if (*p a oe rat) 


Of 


continue; 

strcpy (if _name,p); 

if ((yyin-fopen(if name, "r")) == NULL) { 
fprintf(stderr, "mvl: can’t open $s\n",if name) ; 
exit (1); 

) 

init stats (&Dtt_ stat); 

init stats (&PA_stat); 


/* Set the EOF context */ 
if (setjmp(eof_ context) == 0) ( 


/* While still expressions to parse in this file ... */ 
while (expr() > 0) ( 

nvar = £ orig.nvar; 

radix = £ orig.radix; 

tot_terms t= £ orig.nterm; 

tot exprtt; 


if (D flag) { /* do the Dueck © Hitlers 
- Dueck Miller(); 
FINAL = DM; 
) 
if (F_flag) ( /* do the Pomper & Armstrong */ 
Fomper Armstrong (); 
if (!D_flag) 
FINAL = P_A; 
) 
if (G_ flag) ( /* do the Gold (which chooses between above) 
Gold (); 
if (of_name(0}) /* output the final expression */ 
output expr (HEUR) ; 
FINAL = HEUR; 
) 
else [ 
if (of name({0)) ( 
if (D flag) /* output the final expression */ 
output expr (D_M); 
else if (PF flag) 
output expr (F_A); 
) 
J 
verify (); 
dealloc_expr (&® orig); 
) 


) 
else if (!tcompleted) 
syntax_error (ERR_UNEXF_ EOF); 


/* Print the statistics */ 
if (9 flag) ( 
if (D_flag) 
print DM stats(); 
if (P_ flag) 
print PA_stats(); 
) 
fclose (yyin); 


J 
if (of _name[0]) { 


/* Build the output file for the layout generator by appending 
* the list of final expressions in the temporary file to 
a header consisting of radix, inputs and outputs. 


Final format is: 


* 
& 
* 
* 
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a 


radix inputs(nvar) outputs (tot_expr) 
coeff 

lower upper 

lower upper 


99 (sentinel) 
coeff 

lower upper 
lower upper 


99 (sentinel) 


etc 


> + © © & + FF F&F &F S&F SF SF HS FD 


/ 


/* Close the temporary file, reopen for reading and open 
* the final output. 
dT 

fclose (yyout);? 

if ((yyin = fopen(tmp_name,"r")) == NULL) { 


fprintf(stderr, "mvl: can’t open temp file - %s\n",tmp_name) ; 
exit (1); 

} 

if ((yyout = fopen(of_name,"w")) == NULL) ( 
fprintf(stderr, "mvl: can’t open output file - %s\n",of name); 
exit (1); 


) 


/* Build and output the header */ 
sprintf (expression, "%d %d %d\n", radix, nvar,tot_expr); 
fputs (expression, yyout) ; 


/* Append the temp file, then delete it */ 

while (fgets(expression,MAX BUF,yyin) != NULL) 
fputs (expression, yyout); 

fclose (yyout); 

fclose (yyin) ; 

unlink (tmp_name) ; 


99 


Gold () 


:function: 
- Implement the Gold heuristic (which chooses between F&A and D&M, 
preferring D&M on ties). 
:globals: 
E finalf{) 
F orig 
HEUR 
:called by: 
main () 


char * pe 


if (£ final(0 M).nterm < E final(F Aj}.nterm) [ 
HEUR = D_M; i “ 
p = "D&M"; 
) ; 
else if (FE final(D_M}].nterm > E_ final(P_A].nterm) ( 
HEUR = F A; 
p= "PGA": 
} i] 


else ( 
HEUR = D M; 
Pp = Weems 


) 
1£ “(!q flag) °( 
printf£(" Gold(%s): %4d/%-4d implicants - %$4.2£%%8\n\n", 
p,E— final(HEUP].nterm,£ orig.nterm, 
(float)E final[{HEUR] .nterm/ (float) BE orig.nterm); 


60 


print map () 
a a ee ee a ee ee 
:function: 
- Print the Karnaugh map of E work in its present state 
:globals: 
F work 
:called by: 
Dueck Miller () 
Pomper_ Armstrong () 
:calls: 


register , hag |e 
int X{MAX_VAR+1); 
int *\- 


for (40; £ < nvar; i4+) X{[i}) = 0; 
for (1"0; i < nvar;) ( 
Vea aneval (GE work, Xx); 
printé£ ("%3d%c",V[EVAL),V[HLV)?".%3% "); 
Xfij++; 
for {71° < nvar;) ([ 
4£ (X[i}) >= radix) ( 
X{i}) = 0; 
if (1 < 2) 
print£("\n"); 
L443 
Xfij+4; 
) 
else ( 
4 = O; 
break; 


61 


print _implicant (s,X,1) 


char #9; 

int *Xe 

Term *Is 

| ee 


:function: 
- Print the Most Isolated Minterm X and the implicant selected 
to cover it I. 
scalled by: 
Dueck Miller) 
Pomper Armstrong () 


register i; 


printf(" %s MIM: %2d",3s,X[nvar)); 
for (i1=0: 1 < nvar; i++) 

printé ("*x%d(%$2d, $2d)",i4+1,X[(i),X iJ): 
print £("\n")> 


printf (" Imp: %2d", i->coctket) ; 
for (1=0° i < nvdar; 1144) 
printf ("*xX%d ($2d, $2d) ",i+1, I->B(1) .lower, I->B(i} .upper); 
printf ("\n\n"); 
} 


print terms (£) 
Expression *E; 


sfunction: 

- Print the expression pointer to by £& 
scalled by: 

Dueck Miller () 

Pomper Armstrong () 


into oe 


printf ("\nExpression:\n\n"); 
printf(" radix: %2d\n",E->radix); 
peintf(" mnvars: %2d\n",E->nvar); 
printf(" nterms: %$2d\n",E->nterm); 
for (1-0; i < E->nterm; itt) (¢ 

printf(" coef£: %d\n",E->T[i}).coeff£); 

for (j-0? j < E->nver; j++) 

printf (" Xd (Sd, $d) \n", 
jt1,f£->T(1).8(3).lower, E->T[1i].B(j).upper) ; 

) 
printf ("\n"); 
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output expr (heur) 
int heur; 
eee oe ee = a ee ee 
:function: 
- Append the final expression to the temporary file. 
- heux indexes into E_ final. 
-~ If the final expression is not at least as good as the original 
input expression, output the original. 
- Output format as follows: 


coeff 
lower upper 
lower upper 
99 (sentinel) 
:globals: 
yyout 
stdout 
BE orig 
FE final 
scalled by: 
main () 


register i, 4; 
Expression *E; 


if (yyout t= stdout) { 
if (BE orig.nterm < E_ finalf{heur].ntexrm) 
E= &E oxig; 
else 
Ee GE finalf{heur]; 
for (i=0; i < E->nterm; itt) { 
fprintf (yyout, "d\n", E->T[i].coeff); 
for (j=0; j < E->nvar; j++) 
fprintf (yyout,"%d %d\n", 
E->T[i).B[j].lower, 
E->T[i).B[j}).upper); 
) 
fprinté (yyout, "99\n"); 
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print source (E,3s,n) 
Expression *8; 
char “3s; 
Je ----------2- 2-25-88 
:function: 
- Frint the expression pointed to by E to the file xf_name in 
the form recognizeable by the parser. 
:globals: 
xf name 
:called by: 
Oueck iter) 
Pomper Armstrong () 


reaister i, 4; 
FILE bald 


if ((£ = fFopen(xf name,"a")) ©= NULL) { 
fprintf(stderr, "mvl: can’t open $s\n",xf name); 
return; 


fprint£(£,"# %3: %d/%d\n",3,n,E->nterm) ; 
Eprint£&£(£, "td: %d:", £->radix, E->nvar); 
for (1=<0; i < E->nterm; it+) ( 
forint’ (t, "\n 4%d", E->T[1}.coeff); 
for (j=0; 4} < B-Snvar; 4j++) ¢ 
Forint f (f£, "*xtd ($d, $d)“; 
341, 8->T[(i}) .B[ 4] lower, E->T[iJ.B(4).upper); 
) 
) 
Forintl(l, "2 \n"): 
fclose(f); 


init stats (S) 
MVL stats os 


/* 


) 


i ed 


-function: 

- Initialize statistics structure pointer to by S$ 
:called by: 

main () 


Eotsexpr = 0) 

tot terms = 0; 

S->calls eval = (0; 
S=>calls cl = 0; 

S=-calls mim = 0; 

S->calls gen bounds = 0; 

s- ->calls pick_ implicant = 0; 
Sa calls  nexcnimplicant = 0; 
S->calls _valid _implicant = 0; 
S- >calls_compute_rbc = 0; 


print _PA_stats() ; 


/* 


—— —_ 


:function: 

- Print the stats on the P6A heuristic 
:globals: 

tot terms 

POU ,expr 

radix 

nvar 

PA_stat 
scalled by: 

main () 


— SD ee Ge Gee ee ee ee ee ee ee ee ee ee ee ae a ae oe ee Gee ee ee ee ee ee ee ee oe = ee eo oo oo om oo oo oo © om oe © © oe = om oe om om om oo @ = == 


printf ("Statistics for P&A:\n\n"); 


printf (" Nim Num Num Avg Num\n"); 

printf ("Radix Var Expr Terms Texrm/Expr\n") ; 

printf£("t4d %4d %6d 86d %$8.2£\n\n", 
radix,nvar,tot_expr,tot_terms, (float)tot_terms/(float)tot expr); 


printf (" Eval Fick Gen Next Valid\n"); 
printf (" Expr MIM Bounds Impl Impl \n"); 
print£ ("Tot: S*llild %81d &8ld %@ld %&8ld\n", 


FA _stat.calis eval,FA_stat.calls _mim, 
EA_ stat.calls | agen pounder FA_stat.calls next_implicant, 
FA stat. calls valid_ implicant); 
printf£("Avg/Expr: %11.2£ %$8.2£ %$8.2£ %$8.2£ $8.2£\n", 
(float) FA_stat.calls eval/(float)tot_expr, 
(float)FA_stat.calls mim/(float)tot_expr, 
(float) PA_ stat.calls gen_ bounds/ (float) tot _@Npr, 
(float) PA_ stat.calls _next _implicant/ (float) tot _@xXpr, 
(float)PA_ stat.calls _valid _implicant/ (float) tot expr); 
printf£("Avg/Term: %11.2f %38.2f %$8.2£ $8.2£ %8.2£\n\n", 
(float)PA_stat.calls eval/(float)tot_terms, 
(float)PA_stat.calls _mim/ (float) tot _terms, 
(float) PA_ stat.calls gen_ bounds/ (float)tot _texms, 
(float) PA stat.calls _next _implicant/ (float) tot _terms, 
(float) PA stat. calls _valid _implicant/ (float) tot_ terms); 


65 


print DM stats () 
[®t oncn- onan ccmees See en ae ae mee om ee ee ee ee 
‘function: 
- Frint the stats on the D&M heuristic 
:globals: 
tot terms 
tot _expr 
radix 
nvar 
DM stat 
scalled by: 
main () 


printf ("Statistics for D&M:\n\n"); 


print (" Num Num Num Avg Num\n"); 
printf("PRadix Var Expr Terms Term/Expr\n"); 
printf("%4d %4d %6d 86d %$8.2f\n\n", 
radix,nvar,tot_expr,tot_terms, (float)tot_terms/ (float)tot expr); 
princrl. , Eval Comp Pick Comp Gen Next 
Prince: Expr Cr MIM RBC Bounds Impl 
printf ("Tot: Blild BAld S8Bld #Bld $Ald &Bld SAId\n", 
Dit stat.calls eval,DM stat.calls cf,DM stat.calls mim, 
DM stat. calls compute rbc, 7 - 7 
DM stat.calls gen_bounds,Dt_stat.calls next implicant, 
DM stat. callss valid _implieant) + 
printf ("Ava/Expr: %11.2£ %8.2f 88.2f %8.2F %$8.2f $8.2 %8.2£f\n", 
(float) DM. stat calls eval/ (float) tot (expr, 
(float)DM stat.calls cf/(float)tot expr, 
(float) DM_stat.calis mim/ (float)tot_ expr, 
(float)DH_stat.calls_ compute rbc/(float)tot_expr, 
(float) DM stat.calls gen _bounds/ (float) tot _expr, 
(float) DM stat-calls. next _implicant/ (float) tot _eOxpr, 
(float)DM stat.calls_ valid _implicant/ (float) tot expr); 
printf ("Avg/Term: %11.2f %&8. 2f &8.2£ &8.2£ 88.2£ $8.2£ %8.2£\n\n" i 
(float) DM stat.calls eval/(float)tot_ terms, 
(float) Di _stat.calls cf/ (float)tot terms, 
(float) DM stat.calls mim/ (float)tot terms, 
(float) DM stat.calls compute rbc/(float)tot_terms, 
(float)Dt! stat.calls_ gen_ bounds/ (float) tot_ terms, 
(float) DM stat.calls_ next _implicant/(float)tot _terms, 
(float})DM_stat.calls_ valid _implicant/ (float) tot _terms); 


fatal (s) 

char sal By 

( 
fprintf(stderr, "&s\n", 8); 
exit (1); 
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Valid\n 
Impl\n" 


/* defs 
- E 
a] 


fincliud 
finclud 


"define 
fdefine 
fdefine 
fdefine 


fdefine 
fdefine 


#fdefine 
fdefine 


fdefine 
fidefine 


-h 


xternal definitions for mvl 


e <stdio.h> 
e <set jmp.h> 


MAX BUF 4096 
HAX PATH 64 


MAX VAR 32 
MAX_INT 32765 
YES 1 

HO 0 

EVAL 0 

HLV 1 

DM 0 

P_A 1 


/* parser definitions and data structures 


fdefine 
fdefine 
fdefine 
fdefine 
"define 
fdefine 
fdefine 
fdefine 
fdefine 
fdefine 
fdefine 
fdefine 
#idefine 


fdefine 
fdefine 


fdefine 
fdefine 
fdefine 
fdefine 
fdefine 
fdefine 
fdefine 
fdefine 
fdefine 


extern 
extern 


extern 
/* MVL 
typedef 
typedef 
typedef 
typedef 


struct 


TOK EOF 0 

TOK NUMBER 2 

TOK VAR 3 

TOK COLON 4 

TOK STAR 5 

TOK _LFAREN 6 

TOK COMMA 7 

TOK RPAREN 8 

TOK SIGNED 9 

TOK DONTCARE 10 
TOK SEMI 11 

TOK PLUS 12 

TOK MINUS 13 
TOK_EOL 14 
TOK_SYNERR -1 
ERR_RADIX RNG -2 
ERR_NVAR_RNG -3 
ERR COEFF RNG -4 
ERR_VARINDEX RNG -5 
ERR_DUF_VAR -6 
ERR LOWER RNG = -7 
ERR _UFFER RNG -8 
ERR_UFFER LOWER -9 
ERR_UNEXP EOF -10 


char expression(MAX_BUF+1]; 
int token _type; 


FILE *yyin, *tyyout; 

definitions and data structures */ 
struct expr struct Expression; 
struct term struct Term; 
struct bound struct Bound; 


struct term struct Implicant; 


expr struct [ 
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em mee eee ec cee ce ce ee cee ee ee ee ee ee ee ee ae ee 


int radix, 
nvar, 
nterm; 

Term id IES 


es 
’ 


struct term struct 4 
int coeff; 
Bound *B; 

I; 


struct bound struct [ 
int lower, 
upper; 


be 
’ 


extern Expression 
FE orig, 
E work, 
EF final[2); 


extern int HEUR, FINAL; 


extern int 
nvar, 
radix; 


Perm *alloc term(); 
Bound *alloc_bound(); 


extern char 
if name([ttaxX PATH+1), 
xf name [tix PATHHI], 
Of name [MAX_PATH#+1); 


extern tfiag, 
e flag, 
1 Jf laG; 
x flag, 
m flag, 
G flag, 
v_flag, 
q_ flag; 


extern FILE *xfile; 
extern int completed; 


extern 
Int Cot expe, 
tot terms; 


struct mvl_ stats [ 
long 

calls eval, 
Cali- ce, 
calls mim, 
calls _gen_bounds, 
eslis next implicant, 
calls pick_implicant, 
calls valid implicant, 
calls compute rbc; 


); 
typedef struct mvl_ stats MVL_ stats; 


extern MVL stats 
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*STAT, 

DM stat, 
EA stat, 
G stat; 


extern float FO ratio; 


int *mim(); 

int *next_coord(); 
Bound *gen_bounds(); 
Term *next implicant ()? 
Term *pick implicant (); 
int *eval(),* eval(); 
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/* parse.l 


- lex specification file for the input scanner/parser for the MVL 
parcser/compiler 


‘ 


expression grammar: 


expr i= 
radix ':’ num_var ‘3! term list 


term _list := 
term °3° 
| term ‘'+° term list 


term := 
coeff ‘*" var list 


var list := 
variable 
| variable '*’ var list 


variable := 
var id * (" Timit ' |” timtee” 


verivication test grammar: 


verification :=- 
variable °;’ 
| variable verification 


constraints: 


radix >= 2 

num var >> 1 

1 <= coeff <= radix (radix == don’t care) 

0 <= limit <= radix-1 and lower Limit <= upper limit 
ex? <= var id <= x (nimevary® a = 


data structures 
expression 

radix } 

nvar j 


nterm } 
term * j 


—_— = = ae 


~ An expression structure contains the radix of thea logic 
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td 


expression, the number of variables, the number of variables 
and a pointer to an array of terms 


term 


{ coeff J} 
{ bounds *] 


- A term structure contains a coefficient and a pointer to a 
list of bounds (length equal to the number of variables 
allowed in the expression) 


bounds 


{ lower J[ upper J 
{ lower J[ upper J] 


e 


- A bounds list is an array of pairs (lower and upper bounds) 
which may range over the radix of the expression, one pair 
for each possible variable in the term j 


finclude "“*defs.h" 


*) 
8 


"e"E*\nj*\n ( 


) 


{\t 
} 


an | 


[O-9)+ [ 


strcat (expression, yytext) ; 
return (TOK_NUMBER) ; 


PSI te-3)+ 


) 


apr 


) 


streat (expression, yytext); 
return (TOK_SIGNED) ; 


strcat (expression, yytext) ; 
return (TOK_DONTCARE) ; 


[xX] {1-9} (0-9)* [ 


} 


Han 


strcat (expression, yytext) ; 
return (TOK_VAR) ; 


{ 


strcat (expression, yytext); 
return (TOK_COLON) ; 


( 


Strcat (expression, yytext)?; 
return(TOK_ STAR) ; 


(al 


sak Sacaae | 
streat (expression, yytext); 
return (TOK _ LPAREN) ; 


) 


win ( 
streat (expression, yytext); 
return (TOK_FLUS) ; 


streat (expression, yytext); 
return (TOK MINUS); 


R a { 
streat (expression, yytext); 
return (TOK COMMA); 


"yn" { 
strcat (expression, yytext) ; 
return (TOK_REAREN) ; 


return(NO); 


en { 


streat (expression,yytext); 
return (TOK_SEMI) ; 


f NENEAN\ ET 4 
expression[0) = '\0’; 

} 

$% 


/* parser NT functions for expressions */ 


expr {) 


{ 


int rad,nvar; 


init vexrpe(); 

expression[0} = ’\0’'; 

1f (lyytext[O) Jf match(";")) 
next token(); 

completed = 0; 

if (token type == TOK_FOF) 
return({0); 


if (!radix_val{)) { 
syntax error (TOK_SYNERR); 
return (NO); 

) 

sscanf (yytext, "td", &rad); 

if {rad < 2) { 
syntax error (ERR_RADIX_RNG); 
return (No); 


) 


E orig.radix = rad; 


next token(); 
if ({Imatch(":")) { 


fe 


syntax error (TOK_SY¥NERR); 
return (NO); 
) 


next token(); 

if (tnum_var()) ( 
syntax error (TOK_SYNERR); 
return (NO); 


) 


sscanf (yytext, "$d", &nvar); 

if (nvar < 1) { 
syntax error (ERR_NVAR_RNG); 
return (NO) ; 

) 


FE orig.nvar = nvar; 


next token(); 

if (tmatch(":")) ¢ 
syntax error (TOK_SYNERR) ; 
return (NO)? 

) 


next token(); 

if (ftterm list()) [ 
syntax erro -(TOK_SYNERR) ; 
return(NoO); 

) 


1f (tmatch(";")) { 
syntax error (TOK_SYNERR); 
return (NO); 
) 
completed = 1; 
return(YES); 
) 


radix val() 
( 
1£ (token type == TOK_NUMBER) 
return (YES); 
return(No); 


) 


num_var () 
( 
if (token type == TOK_NUMBER) 
return(YES); 
return (NO); 
) 


term list {) 
( 
register i; 
1£ (tterm()) 
return (NO); 
1£ (match("4+") | match("-")) ¢ 
next token (); 
while (term_list()); 
) 
else 
while (term list()); 
for (1-0; i < E_orig.nvar; 144+) 
14£ (6 orig.T(B& orig.nterm-1].B(i]J.lower == -1) 
FE orig.T(£_orig.nterm-1}.Bli].lower = 0; 
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return(YES); 
) 


term () 
( 
int coeff; 
if ((token_type == TOR_NUMBER) || (token_type == TOK _SIGNED)) ( 
sscanf (yytext, "%d", &coeff); 
VE (tie oi Lag) it 
if ((coeff < 1) || (coeff >= E_orig.radix)) ( 
syntax error (ERR_COEFF RNG); 
return (NO) ; 


) 


else if (token_type == TOK_DONTCARE) ( 


coeff = E_ orig.radix; 
) 
else ({ 

return(NO)?; 
) 


E orig.ntermtt; 
EZ orig.T = alloc term(E_orig.T, coeff,E orig.nterm); 


next tokén(); 

1£ (!ftmatch("*")) [ 
Syntax error (TOK _SYNERR) ; 
return{(NO); 

} 

next token({); 

1£ (!var Fist {)) 
return (NO); 

return(YES); 


var list () 


if (!variable ()) 
return (NO); 
next token (); 
if (mateh(**%)) *{ 
next token (); 
while (var list ()); 
) 
return(YES); 
) 


variable {) 
( 
int varindex, lower, upper; 
if (!var_id()) 
return (NO); 
sscanf (&yytext[1}), "%d", &varindex); 
4f ((varindex < 1) || (varindex > E_orig.nvar)) ( 
syntax error (ERR_VARINDEX_ RNG); 
return(NO); 
) 
varindex--; 
if (E_orig.T{E_orig.nterm-1).B[varindex].lower {= -1) ( 
syntax error (ERR_DUP_VAR);? 
return (NO); 
! 


next token(); 
if (tImatch("(")) ( 
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syntax error (TOK_SYNERR) ; 
return (NO) ; 

) 

next token(); 

4£ (fLimdit()) ( 
syntax error (TOK_SYNERR); 
return (NO); 

) 

sscanf (yytext, "%d", &lower); 

1£ ((lower < 0) ){ (lower > E_orig.radix-1l)) { 
syntax error (ERR_LOWER_RNG) ; 
return (NO); 

) 


next token(); 

1£ (!tmatch(",")) ( 
syntax error (TOK_SYNERR); 
return (NO); 

) 

next token(); 

4£ (tLimit()) { 
syntax error (TOK_SYNERR); 
return (NO); 

) 

sscanf (yytext, "%&d", &upper); 

1£ ((upper < 0) ) (upper > E_orig.radix-1)) ( 
syntax error (ERR_UFFER_RNG); 
return (NO); 

) 

1£ (upper < lower) ( 
Syntax error (ERR_UFPPER_LOWER) ; 
return (NO) ; 


) 


next token(); 

1£ (tmatch(")")) f{ 
syntax error (TOK_SYNERR) ; 
return(NO); 

) 


EF orig.T({E orig.nterm-1}].B[varindex}].lower = lower; 
E_orig.T[E_orig.nterm-1].B(varindex].upper = upper; 
return (YES); 

) 


var id() 


( 
if (token_type == TOK_VAR) 


return (YES); 
return (NO); 
) 


limit () 


( 
if (token _type == TOK_NUMBER) 
return(YES); 
return (NO); 
) 


/* verification NT functions */ 


verify () 
( 


register {; 
Bound B([MAX_VAR+41); 
Term a 
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int first, first time = 1; 
int “value, *X; 


for (77) ( 
1£ (match (";")) 
next token(); 
T.B = B; 
for {(i1=0O; i < nvar; itt) 
B{ij.lower = -1}; 
if (Iver variable (B)) 
return; 
completed = 0; 
while (ver variable (B)); 
4f (!Imatch(";")) 
syntax error (TOK_SYNERR) ; 
completed = 1; 
if (!v_ flag) 
continue; 


first = 1; 
4f (first time) ( 
printf ("Verification:\n\n"); 
first time = 0; 
) 
for (i-0; i < nvar; i144) 
printf£(" x%d(%d, $d) \n",41+1,B(i). lower, B[i] .upper); 
Printetsvnc 
while ((X = next_coord (ST, first)) != NULL) ( 
first = 0; 
value =~ eval (6&E final [(FINALJ,X); 
printf(™  (*); 
for (i=0: i < nvar; itt) ( 
preintt (*td", x{(i)0-> 
4¢ (i41 < nvar) print£t£(", °): 
) 
printf£(") = %3d\n",value[EVAL])); 


] 
printf ("\n"); 


) 


ver, Variab1e:(8) 
Bound “Be 
{ 
int varindex, lower, upper; 
if (lvar_id()) 
return(NO); ; 
sscanf (&yytext (1), "td", &évarindex) ;: 
if ((varindex < 1) [|] (varindex > nvar)) { 
Syntax errer(ERR_VARINDEX RNG); 
return (HO) ; 
} 
varindex--; 
if (Bivarindex].lower != -1) [ 
Syntax error (ERR_DUP_VAR); 
return(NoO) ?; 
) 


next token (); 

4£ (Imatch("(")) ¢ 
syntax error (TOK SYNERR) ; 
return(NoO); 

) 

next token(); 

Lf (flimit()) f{ 
syntax error (TOK_SYNERR) ; 


return (NO) ; 
} 


ascanf (yytext, "%d", Slower); 


if ((lower < 0) {| (lower > radix-1)) ( 
syntax error (ERR_LOWER_RNG); 


return (NO); 


) 


next token(); 
1£ (tmatch(",")) ( 


syntax error (TOK_SYNERR); 


return(NO); 
) 
next token(); 
1£ (tLlimit()) ¢ 


syntax error (TOK_SYNERR); 


return {No) ; 
} 


sscanf (yytext, "%d", supper) ; 


1f ((upper < 0) [|] (upper > radix-1)) ( 
syntax error (ERR_UFFER_RNG); 


return(NO); 
) Z 
if (upper < lower) ( 


syntax error (ERR_UFPER_LOWER) ; 


return(NoO); : 


) 


next token(); 
if (fmatch(")")) ( 


syntax error (TOK_SYNERR) ; 


return(NO); 
} 
B(varindex] .lower = lower; 
B(varindex]) .upper = upper; 
next token(); 
return (YES); 
) 


/* scanner utility functions -- 


syntax error (code) 


( 


static char *errors[] = ( 


/* TOK_SYNERR dy 
/* ERR_RADIX_RNG acl 
/* ERR_NVAR_RNG 77 
/* ERR_COEFF_RNG ai 
/* ERR_VARINDEX RNG */ 
/* ERR_DUP_VAR */ 
/* ERR_LOWER_RNG as 
/* ERR_UFPER_RNG Ar 
/* ERR_UPFER_LOWER */ 
/* ERR_UNEXP_EOF ei 


? 


code = abs(code); 


"mH 
¢ 


"syntax error", 


“radix out of range (>= 2)", 

"# of var out of range (>= 1)", 

"coeff out of range (1 <= coeff <= radix-1)", 
"var index out of range (1 to # of var)", 
"duplicate variable id", 

"lower bound out of range (0 to radix-1)", 
“upper bound out of range (0 to radix-1)", 
“upper bound must be >= lower bound", 
"unexpected end of file", 

NULL 


printf ("%s\nts\n",errors(code],expression); 


exit (1); 
) 


next_token () 


(tl 


token type = yylex(); 

if (token_type == TOK_SYNERR) { 
Syntax_error (TOK_SYNERR) ; 

) 

return (token type); 


) 


match (9) 
char .s; 
( 
Z£ (!stremp(s,yytext)) 
return (YES); 
else 
return (NO); 
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/* dm.c 


- This module implements the Dueck & Miller heuristic for minimizing 
an expression. 


s/ 
finclude "defs.h" 


Dueck Miller () 
/* ae ee ee ee ee le ee ea a oa oe eee 
¢:function: 
- Ferform the Dueck & Miller heuristic on the input expression 
salgorithm: 
Start with a working copy E_work of the original function E_ orig; 
Initialize a final function Ff final; 
While (there are still minterms to pick) ( 
Pick a minterm X from E_work; 
Pick the best implicant I for xX; 
Subtract I from E_ work; 
Add I to E_ final; 
) 
:globals: 
E orig 
e flag 
m flag 
q flag 
G_flag 
FO ratio 
sside effects: 
STAT 
HEUR 
E work 
E final(} 
:called by: 
main () 
:calls: 
dealloc expr () 
dup expr () 
print terms (} 
print map (} 
mim (} 
pick implicant () 
subtract implicant () 
print source () 


——-—ae eS Se SS Se eS eS ee 2 Se ee ee eS ee Se ew ee ew 2 2 Se Se Se SP SSO SS OS SF SF SF SF SO SF SS SSF SF SF SF SSO FSF OOF e So we wH— ww See SH SH oe 


int num_imp1 = 0; 
int *x; 
Term ai 
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float ratio; 


if (B_final{D_M}.T f= NULL) 
dealloc_expr(GE_ final{D_M)); 


STAT = &DM stat; 

HEUR = DM; 

dup_expr(&B_work, GE orig); 

E final({HEUR].nterm = 0; 

E final {NEUR}). radix ~ £ orig.radix; 
E final{NEvuR}.nvar = E orig.nvar; 

E final{HEUR).T = NULL; 


if (e flag) 
print terms (&E orig); 

if (m_ flag) ( a 
printf£(" Orig map (D&M):\n"); 
print _map(); 


£Or Aes) 
2 4f ((X = mim(&E work)) == NULL) 
break; 
I = pick implicant (X); 
num_impl+#+; 
subtract implicant (I); 
if (4 flag) 
print implicant ("D&M",X,I); 
if (m_flag) 
print _map(); 
) 
ratio = ((float)num_imp]/ (float) E_orig.nterm); 
ff (FOUrsttorl=20-0)42| 
if (ratio > FO_ratio) { 
print _source (&E orig, "D&M",num_imp]); 
} 
} 
if (!q flag && !G flag) ( 
printf(" D&M: %4d/%-4d implicants - %4.2£%%\n\n", 
num_impl,& orig.nterm, ratio); 
) 
dealloc expr(&E work); 
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static int cf (E,x) 
Expression *E; 
int x; 


static int _cf£(E,X) 


:function: 
- Determine the clustering factor for (x1,x2,x3,x4....xMm) 


- Builds a vector of coordinates for cf () 
- Local to dm.c 
:calls;: 
_cf () 
sreturns: 
~ An integer clustering factor 


int *X; 


X = (int *) &x; 
return( cf (E,X)); 


Expression *£; 
int *X; 


/* 


:function: 
- Compute the clustering factor at X, where X is a vector of 
coordinates. 
- Local to dm.c 
:globals: 
radix 
nvar 
:side effects: 
STAT 
scalled by: 
mim () 
:calls: 
_eval () 
veopy () 
sreturns: 
~ An integer clustering factor 


int nterm = F->nterm: 
register i,j,k; 
int value[2], 
expanded, 
vall{2), 
val2(2); 
int cf = QO, 
dea = 0, 
ea = 0; 


STAT=>calls crt; 

vcopy (value, eval (E,X)) 

if ((value[FVAL] < 1) | 
return (MAX_INT); 


! (value[EVAL] >= radix)) 


/* for each variable (direction)... */ 
for (i-0; i < nvar; i4+ ) f 
expanded = 0; 
/* If not on a left hand edge, move left */ 
if (xX(i}) ( 
X[i}--; 
veopy(vall, eval (E,X)); 


8] 


4£ (vall[ EVAL) >= value[EVAL]) ({ 
expanded+#; 
eatt+; 
) 
X[i)++; 
} 
/* if we didn’t start ona right hand edge, move right */ 


Lf (X(i] >< (redix-1))) | 

X[i}+4? 

vcopy (val2, eval (E,X)); 

L£ (val2[EVAL] >= value[EVAL]) | 
expandedt#; 
eat+; 

) 

Se er 


) 
1£ (expanded) 


deatt; 
/* compute the clustering factor */ 
cf = (dea * (radix-1)) + ea; 
return(cf£); 


static int *mim(E) 
Expression *E; 


:function: 
- Find the Most Isolated Minterm in the expression pointed to 


by FE, and return its coordinates as a vector. 
- Local to dm.c 
:globals: 
nvar 
sside effects: 
STAT 
:called by: 
Dueck Miller () 
:calls: 
next coord () 
_eval () 
veopy (} 
_cf () 


sreturns: 
- A vector of integers representing the coordinate of the most 


isolated minterm, or NULL if no more minterms. 
- The value at that location is also returned as the last integer 


in the vector. 


int cur_val = E->radix, 
cur CF = MAX_INT, 
Ge, 
value[2]), 
term; 
int *X,*next coord(),4; 
static int save_coord(MAX_VAR+t1}; 


STAT->calls mimt+; 
for (term=-0; term < E->nterm; termif+) = ( 
i= 1; 
while ((X=next coord (&(E->T[term)),i)}) != NULL) [ 
veopy (value, eval (E,X)); 
if (value[EVAL} > 0) { 
if (value(EBEVAL) < cur_val) [ 
cur val = value[EVAL); 
for (i=0; i < BE->nvar; i++) save coord(i) = X[(iJ; 
Cire Cr = et {5 Xx); 7 
) 
else if (value[EVAL] == cur _val) ( 
c£ = cf£(E,X); 
LE MCE <cur cr) { 
cur CF = cf? 
for (i=—0; i < E->nvar; i++) save coord(i} = xX[(i); 


i = QO; 
} 


) 
if (cur_CF == MAX_INT) 
return (NULL); 
save coord(£->nvar) = cur val; 
return (save coord); 
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static int valid implicant (I) 
Term ban 
[8 2------ 2 on oo ee ee == 
function: 
- Decide upon the validity of implicant 1 
- Local to dm.c 
sglobals:; 
E work 
E orig 
:side effects: 
STAT 
:called by: 
pick_implicant () 
scalls: 
next coord () 
_eval () 
veopy () 
sreturns: 
1 if 8 valid implicant 
0 if not 


5 Mig ofa at 

int init = 1; 

int value = I->coeff£; 
int Vo[2],Vw[2]J; 


STAT->calls valid implicant+#; 
while ((X = next _coord(I,init)) [= NULL) { 
init = 0; 
veopy (Vw, eval (&E work,X)); 
veopy (Vo, eval (&E orig,X)); 
Lf (((Vw[EVAL] < value) && IVw[HLV]) && (Vo[EVAL] < (radix-1))) 


return(0); 


) 


return(1); 
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static int compute_rbc(I) 
Term al G2 


:function: 
- Compute the RBC for the given implicant 
= Local fo dm.c 
:globals: 
radix 
nvar 
7sidevetteccs: 
STAT 
scalled_by: 
pick_implicant () 
seals: 
next coord () 
_eval() 
veopy ({) 
sreturns: 
- an integer RBC 


int *X; 
int I value = I->coeff}; 
register 1; 
int value[2], 
neighbor value[2], 
good, 
bad, 
rbc = 0, 
init = 1; 


STAT->calls_ compute _rbctt; 
/* for each coordinate in the implicant ... */ 
while ((X = next _coord(I,init)) != NULL) [ 
init = 0; 
veopy (value, eval (GE _work,X)); 
1£ (value[{EVAL) == radix) 
continue; 
/* for each direction ... */ 
for (i-0; 1 < nvar; i++) ( 
good = Q; 
bad = Q; 
if (value[EVAL) == I_ value) 
good = 1; 


/* if there is a left neighbor, examine it */ 
Lf ((Xfi} t= 0) && (XL) == I->Bli).lower)) { 
X[i)-=; 
veopy (neighbor value, eval (&6£ work, X)); 
X[ i) ++; 
1£ (neighbor value[EVAL]) != 0) ( 
if (neighbor value[EVAL) == (value[EVAL]) - I_value)) 
good = 1; 
if (neighbor value[EVAL] == value[EVAL)) 
bad = 1; 


) 


/* if there is a right neighbor, examine it */ 
1£ ((X{1) t= (radix-1)) 6&6 (x[{1J ==— I->B{i].upper)) [ 
XL} +4+; 
vcopy (neighbor value, eval (&E_ work,X)); 
X{i}--; 
1£ (neighbor value[EVAL] != 0) ( 
if (neighbor value[EVAL] == (value[EVAL] - I_value)) 


4 
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good = 1; : 
if (neighbor value[EVAL] == 
bad = 1; _— 
} 


) 
/* update the rbc */ 


rbc = (rbc - good) + bad; 
) 
) 


return (rbc); 


static Term *pick implicant (X) 
int ie 


:function: 
- Fick the best implicant for minterm X 
:globals: 
radix 
:side effects: 
STAT 
:called by: 
Dueck Miller () 
scalls: 
init implicant () 
gen_bounds () 
next implicant () 
_eval() 
vcopy () 
compue_rbc() 
copy implicant () 
valid implicant () 


sreturns: 
-~ A pointer to a term representing the best implicant. . 


int cur _rbc = MAX_INT, 
rbc = 0; 
Implicant eT: 
static Bound I_bound[(MAX_VAR+1]; 
static Term I_best; 
Bound *B; 
int V{2}; 


STAT->calls_ pick _implicant++; 
I_best.B = I bound; 
init implicant (X); 
B = gen_bounds (X); 
while ((I = next implicant (B)) f= NULL) ( 
veopy(V, eval(&E orig,X)}: 
if (V[EVAL}] == (radix-1)) ( 
for (I->coeff = X[nvar]; I->coefE < radix; (I->coeff)+t+) (¢ 
if (valid implicant (I)) ( 
rbe = compute_rbc(I); 
Pf (xbe < cur rhe) { 
Curerbe = rhc; 
copy_implicant (6&1 best,1); 


} 
) 
else | 
I->coeff = X[nvar}; 
if (valid implicant (I)) { 
rbe = compute rbc(I); 
1fE (xbe < cur rbej 
Cur rbc = rhe; 
copy_implicant (&I_ best,1); 


) 
} 
return(&I_ best}; 
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/* pa.c 


- This module implements the Fomper-Armstrong heuristic for minimizing 
an expression. 


ad 
finclude “defs.h" 
Fomper Armstrong () 


:function: 
- Perform the Fomper & Armstrong heuristic on the input expression 
:globals: 
Fo_ratio 
E orig 
e flag 
m flag 
q_ flag 
4d flag 
G flag 
:side effects: 
STAT 
HEUR 
BE work 
BE final() 
scalled by: 
main () 
scalls: 
dealloc expr () 
dup_expr () 
print terms () 
print map () 
mim () 
pick implicant () 
subtract implicant () 
print source () 


ae eee ee ee ee ee es ee es es ee ee es es es se es es ee es es es es ee es es ee ew es es es a es ee es a es es es ee ee ee ee we we we we es ee ee ee we es ee ee ee eo 


int num_impl = 0; 
int *X:3 

Term a 

float ratio; 


4f (& finalf{P_A}.T te NULL) 
dealloc expr(&E_final[P_A}); 


STAT = GFA stat; 


HRUR = BA; 
dup _expr(&E work,&E orig); 
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E final(HEUR).nterm = 0; 

EB final(HEUR).radix = E_orig.radix; 
E final[(HEUR}).nvar = F orig.nvar; 
BE final{HEUR) .T = NULL; 


if (e_flag) 
print terms (&E_ orig); 

if (m flag) ( 
printf(" Orig map (F&A) :\n"); 
print map (); 


form(s49.{ 
if ((X = mim(&E work)}) == NULL) 
break; 
I = pick implicant (X); 


num impl++; 
subtract implicant (I); 
if (i flag) 
print implicant ("P&A",X,1I); 
if (m flag) 
print map(); 
) 
ratio = ((float)num_imp1/ (float) E orig.nterm); 
if (FOvratio I= 0.0) { 
if (ratio > FO ratio) { 
print source (&E orig, "P&A",num imp]}; 
} 
} 
if (!q flag && !G_ flag) [{ 
printf(" P&A: %$4d/%-4d implicants - $4,.2£%%\n\n", 
num_impl,E orig.nterm, ratio); 
) 


dealloc_ expr (&E_ work); 
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static int *mim(E) 
Expression *E; 


:function: 
- Return a vector corresponding to the coordinates of the current 
Host Isolated Minterm 
:globals: 
radix 
nvar 
:side veffects : 
STAT 
:called by: 
Fomper Armstrong () 
:calls: 
next coord () 
_eval () 
veopy () 
random ({) 
sreturns: 
- A vector of integers representing the current Most Isolated 
Minterm, or NULL if no more minterms. 
- The last integer in the vector is the value of the function at 
that coordinate. 


int weight = MAX_INT; 

int tmp; 

int ex, 41, *nhext coord); 

int first; 

static int save coord(MAX VARt1] = ( —-1e); 
int value[2]; 

int term; 


STAT->calls mimt+; 
for (term = 0; term < E->nterm; termt++) ( 
first = 1; 
while ((X=next coord (&(E->T{term]),first)) t= NULL) ( 
vcopy (value, eval (E,X)); 
1£f (({valuefEVAL]) > 0) && (value[EVAL]) t= radix)) ( 
tmp = random{); 
if (tmp < weight) [ 
for (1i=0; i < nvar;: i++) 
save coord{i] = xX{iJ; 
save _coord([nvar]) = value[EVAL]; 
weight = tmp; 
) 
) 
first = 0; 
) 
) 
1f (weight !e= MAX INT) 
return(save_ coord); 
return (NULL); 
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static int random() 
[Rt -------------- - - - - = = = 5 5 5 5 5 5 5 ee 5 5 eee eee 
S£unceson: 
- Return a number between 0 and 32766 using a linear congruential 
method 
:called by: 
mim () 


me ce ee we we we we we we we ee ee es ss ss es ee es 9 9 ee ee es es ee es ss es 9 9 ee ee we te ee we ee oe oe oe oe oe we ee 


static long a=100001; 


a = (a*125) % 2796203; 
return( (int) (( (float) a/2796203) *32767)); 
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static int valid implicant (1) 
Term hs OF 
[* ann-2-- 2 eo nn oo on n= 
:function: 
- Decide upon the validity of implicant I 
:globals: 
radix 
sigside effects: 
STAT 
:called by: 
pick_implicant () 
tcalls: 
next coord () 
_eval () 
vcopy () 
sreturns: 
num zeroed if valid 
0 if not 


— ee ee ee ee ee ee ee ee ee ee ee ee ee = oe ee = om oe oe 6 oe om om om om oe ce ee ee ee = ee oe oe oe oe = om 6 we = oe oe = oo om om om oe ce Fe ee ee ee ee ee ee ee oe oe 


int *X: 

int init = 1; 

int num zeroed = 0; 
int value = I->coeff; 

int Vo{2),Vw{2); 


STAT->calls valid implicant+#; 
while ((X = next coord (I,init)) {= NULL) ( 
init = 0; 
vcopy (Vw, eval (&E_ work,X))? 
vcopy (Vo, eval (&E_orig,X)); 
1f (((Vw{EVAL) < value) && IVw{HLV]) 6&& (Vo[EVAL] < (radix-l))) 
return (0); 
if (Vw[EVAL] <= value) 
num zeroedt+; 


) 


return(num_ zeroed) ; 
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static Term *pick implicant (X) 
int *X > a 
[®t wonn--~- +--+ ee + ee + en nn no nn oe --- - -- -  e-- 
:function: 
- Pick the best implicant for minterm X 
:globals: 
nvar 
radix 
:side effects: 
STAT 
:called by: 
Pomper Armstrong () 
:calls; 
init implicant () 
gen_bounds () 
next implicant () 
_eval() 
veopy {) 
copy implicant () 
sreturns: 
- A pointer to a Term representing the best implicant. 


int cur _ num zeroed = -1l, 
cur _num_covered ©= -1; 
int num_covered, 
num_ zeroed; 
Implicant sede 
static Bound I_bound{MAX_VAR+t1); 
static Term 1 _ best; 
Pound sa 5 
int vV{2); 
register i; 


STAT->calls pick implicanttt; 
I_best.B = I_bound; 
init implicant (xX); 
B = gen_bounds (X); 
while ((I = next_implicant (B)) != NULL) ( 
veopy{(V, eval (&E orig, X)); 
if (V{EVAL) ==> {(radix-1)) ( 
for (I->coeff = X{nvar); I->coeff < radix; [{I->coeff)++) { 
if (num_zeroed = valid implicant (1I)) { 
num covered = 1; 
for (4-07 i < nvar: i++) 
num covered *@ (1->B[i].upper - I->B[i]).lower + 1); 
if ((num_zeroed > cur _num zeroed) [I ( 
(num zeroed == cur_num zeroed) && 
(num_covered = cur _num_ covered) 
y) { 
cur num zeroed = num_zeroed; 
cur num covered = num covered; 
copy_implicant (&I best,I); 


) 
) 
else { 
I->coeff = X{nvar); 
if (nmum_zeroed = valid implicant (I)) [ 
num covered = 1; 
for (1©0; £ < nvar; i++) 
num covered *= (I->B{i].upper - I->B{i).lower + 1); 
if ({num_zeroed > cur _num_ zeroed) If ( 
{num_ zeroed ©= cur num zeroed) && 
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(num_covered = cur_num_covered) 

ee | 
cur num zeroed = num_zeroed; 

Cur num covered = num_ covered; 
copy implicant (6I_best,I1); 


J 
return(SI_ best); 
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/* comnon.c 


wt 


- Heuristic support functions common to D&M and EEA. 


finclude "“defs.h" 


int *eval (EF, x) 
Fxpression *E; 
ante xe 


/* 


) 


tfunction: 


eee ww we = oe © oe oe eo oe oo oe oe oe © © oe oe oe oe oe © om Oe ee oe ee oe oe ee oe ee ee ee ee ee es es es ee ee es ee et es es es es es 9 


- Evaluate the expression at (x1l,x2,x3,x4,x5....%"n) 
- Builds a vector of coordinates for eval{) 


tcalls: 
_eval() 
freturns: 


- A vector containing the integer result 0 <= eval() <= radix 
and a flag set if the value reached the highest logic value 


int *X; 


X=— (int *) &x; 
return( eval (E,X)); 


int * eval (E,X) 
Expression *£&; 


int 


/* 


xX: 


:£unction: 
- Fvaluate the expression at X, 
:globals: 
nvar 
radix 
:side effects: 
STAT 
called by: 
mim() - pa.c 
valid implicant () - pa.c 
pick implicant() - pa.c 
mim() - dm.c 
valid implicant() - dm.c 
pick implicant () - dm.c 
_c£() 
compute rbc () 
eval () 
gen_bounds () 
print map () 
returns: | 
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where X is a vector of coordinates 


- A vector with the value of the expression at the specified 
coordinate as its first element, and a flag set if this value 
has attained the highest logic value (HLV) 


int nterm = E->nterm; 
register i,j,k; 

int out_of bounds; 

static int v{2]J; 

register rmi = radix-1; 


STAT->calls eval++; 
V{FEVAL) = 0; 
V{HLV) = 0; 
/* for each term... */ 
for (1-0; 1 < nterm; if+) { 
/* €or each variable ... */ 
for (j-0,out_of bounds<0; j < nvar; jtt) ( 
4£ ¢ 
(X(3) < E->T{i).B(j)-lower) |] 
(X{j) > E->T(1i).B({j).upper) 
rot 
out of bounds = 1; 
break? 
J 


J 
1£ (out_of bounds) 
continue; 


/* if this is a don’t care, return the radix */ 
if (E=>T{1].coefl == radix) { 

V{EVAL] = radix; 

return(V); 


) 


V[EVAL] t= E->T{i].coefé£; 

if (V{EVAL}) >= rml1) [ 
/* set a flag which means £ orig was saturated at this X */ 
V{THLV) = 1; 

J 

TE (V{lEVAL)] > cmb) j{ 
V{EVAL] = rml; 

) 

else if (V[HLV] && (V{EVAL] <= 0)) { 
V{EVAL] = radix; 
return (Vv); 

J 

J 


return (V); 
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int 
Term 
int 


*next coord (T, first) 
Ts 
first; 


I I I ee ee ee ee ee ee ee es ll 


:function: 
- Compute the next possible coordinate for term *T 


- If first == 1, initialize the coord vector 


:called by: 
mim() - pa.c 
valid _implicant() - pa.c 
mim() - dm.c 


valid implicant() - dm.c 
compute_rbc () 
returns: 
- An integer vector containing the coordinates. 


static int coord {MAX VAR#1); 
static i? 


/* 4£€ the first time through, load the vector */ 
if (first) [{ 
for (i=-O; i < nvar; itt) [{ 
coord{i} = T->B{i}.lower; 
) 
) 
else {| 
i= 0; 
coord (i) ++; 
LOLs). A 
if (coord[i]) > T->B[i}].upper) ( 
coord{i] = T->B[i].lower; 
it4?: 
if (i >= nvar) 
return (NULL); 
coord{i]++; 
} 
else { 
break; 
J 
J 
) 


return(coord); 


Pound *gen_bounds (X) 
int: 


:funct ton: 
- Generate the permissible bounds around location X in the 
working expression 
:globals: 
radix 
nvar 
E work 
E orig 
:si0e effects: 
STAT 
:called by: 
pick implicant {) - pa.c 
pick_implicant() - dm.c 
[calis: 
evel () 
vcopy () 
s returns: 
- A bounds array 


ee eee ee ee ee oe oe oe om om om om om om om am we ee = Oe Oe ee ee ee we we oe ee ee ee ee ee ee oe om om © om @® oe om be Oe we oe oe we ee we we ee we we we we ie ie ee 


etatic Bound B(MAX_VAR+1]; 
int nterm = F work.nterm; 
register 174, kh: 

int value,Vw{2),Vof{2]}; 

int Xp(MAX_VAR+1}; 


value = X[nvar]; 


STAT=>callsegen boundatr,: 
/* for each variable (direction)... */ 
for (1=0; i < nvar;s 144.) 
/* dup the coordinate */ 
for \[j=02" y= vars tt) Ol) = sla; 
Bli}.lower = X[{1}; 
/* while not ona left hand edge, move left */ 
while (Xpf{i} > 0) [ 
Xp[i)--; 
veoopy (Vw, eval (6E work, Xp)); 
vcopy (Vo, eval (&6E_ orig, Xp))? 
/* Af can’t expand tc lett"... =7 
1f (!( (value > Vw[EVAL]) && (VO[EVAL] < (radix-1)))) { 
B(i} .lower = Xp{iJ; 
) 
else 
break; 
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/* dup the coordinate */ 
for (j=0; 3 < nvar; j++) Xp[3) = X35); 
B(iJ.upper = xX[iJ; 
/* while not on a right hand edge, move right */ 
while (Xp[iJ < (radix-1)) { 
Xp(ijtt+; 
vcopy (Vw, eval (&E work, Xp)); 
veopy (Vo, eval (&E orig, Xp)); 
/* “ff Tecan’ t expand to fight ... */ 
if (!((value > Vw[EVAL]) && (Vo[EVAL] < (radix-1)))) { 
B(iJ.upper = Xp{iJ; 
) 
else 
break; 


) 


return (B); 
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/* Working structures for picking the next implicant within bounds */ 


static Bound IB[MAX _VARt+1}; /* Current bounds */ 
static Term I; /* AImplicant */ 
Staticsint 

I_var, 

I fipet, 

I val; 
int X_orig[MAX_VAR+t1}; /* Where we start */ 


init implicant (X) 
int ie 
[Wwe eww nw enn nn ee nn ee eo ee ee eee ee ee 
:function: 
- Initialize the static term structure above from which successive 
implicants will be returned 
- X is the starting minterm 
sside effects: 
- The structures above 
:called by: 
pick_implicant() - pa.c 
pick implicant () - dm.c 


int nterm & Fy work.nterm; 
register i; 


/* initialize the implicant */ 

I.B8 = IB; 

I.coeff = X[nvar]; 

for (1=0; 1 < nvar: i++) { 
I.B[i).upper = x[i]: 
I.B{i] .lower = X[iJ}; 

] 

I_ var = Qs 

Pificst = 1; 

I val = X(nvar}; 

for (1-0; 1 < nvar? itt) Xcorig(t}) = x [1]- 


Term *next_ implicant (B) 
Bound *R: 
| A ee ee a 
‘function: 
- On each call, return the next implicant within bounds B 
;eide effects: 
STAT 
:called_by: 
pick implicant() - pa.c 
pick _implicant() - dm.c 
returns: 
- An implicant as a term structure 


int nterm = E_work.nterm; 
int Xp (MAX _VAR#1]; 


STAT->calls next_implicant++; 
$f (lo elese) 

I firsc = 0; 

return(&I); 


) 


while (I_var < nvar) { 
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/* expand left */ 
I.B[I_var}.lower-~; 


J/* A£ we can’t go further left, then ... */ 
if (1.B[I_var]).lower < B[I_var].lower) |[ 


/* move back and go right */ 
I.B(I_varjJ.lower = X orig[1_ var); 
I.B(I_var).uppertt; 

/* if we can’t go further right, then ... */ 
if (1.B[I_var]).upper > B[I_var]).upper) [ 


/* reset and go to the next higher dimension */ 
I.B({I var).upper = X orig{I var); 
I_vartt; a ~ 
continue; 
) 
) 


Y var = 0; 
return(&I); 


} 
return(NULL); 
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int copy_implicant (dest, src) 


Term 


— =m Ge Ge Ge Ge Ge Ge @ Ge Ge om oe om om © am ee oe oe Ge om © oe om oe om om om om om om om om om om om om 6 oe oe oe oe ae oe oe oe ae Se ee ee ee oe oe oe oe oe et Ge ee oe oe ae oe = oe 
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*dest,*src; 


:function: 

- Copy the implicant pointed to by src to dest 
:called by: 

pick implicant () - pa.c 

pick_implicant() - dm.c 


register i; 


dest->coeff = src->coeff; 

for {i= 0; i < nvar; itt) [| 
dest->B[i]J.lower = srec->B[i].lower; 
dest->B[(i]J.upper = src->B[{i].upper; 


subtract implicant (1) 


Term 


me ce ee ee ee ee ee ee ee ee ee ee ee ee ee ee ey ee we we ee ee oe Se ee ee ee em oe ee ew ew ee ee ew ee we 
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J 


aT 


:function: 
- Add implicant I to the working expression as a negative term 
(negated coefficient) 
- Add implicant I to tthe final expression 
:globals: 
HEUR 
nvar 
sside effects: 
E work 
FE final{() 


register i,term; 


term = F work.nterm; 
E work.ntermt+; 
E work.T = alloc _term(F work.T, -[(I->coeff),E_work.nterm) ; 
for (i=0° i < nvear: i4t) | 
E work.T[term] .B[i}).lower © I->B[1i].lower; 
E work.T[term) .B{i] .upper = I->B(i] .upper; 
) 


term = — final (HEUR) .nterm; 
E final [("tEUR] .ntermt++3 


& final{HEUR] .T ~ alloc term(E final[(NEUR].T, I->coeff,E final ([HEUR] .ntexrm) ; 


for (1903) 1. = nvars itt) 
E final ("tfUR).T(term] .B(iJ].lower = I->B[{1i].lower; 
E final[(HEUR]).T(term]).B[i]).upper = I->B[i].upper; 


/* veopy |) 


wi 


- copies the value vector from s to d 


veopy (d,s) 
int *d,*s; 


{ 


d(0} = s{0j; 
d{i) = sf{ij; 
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finclude "defs.h" 


/* memory allocation functions 


Term *alloc_term(p, coeff,n) 
Term *p; 
int coeff,n; 


) 


‘function: 


- Allocate space for a term array, initializing the last element 


- If p is NULL, allocate new space 
- If p is not, realloc 

sreturns: 
- A pointer to the Term 


char *malloc(), *realloc(): 
Bound “alloc_bound(); 


if (p == NULL) ( 
1f ((pe(Term *)malloc(sizeof (Term) *n)) "= NULL) 
fatal ("alloc _term(): out of memory \n") ; 
p->coeff — corff; 
p->B = alloc bound); 
} 
else [ 
1f ((p=(Term *)realloc(p, sizeo£ (Term) *n)) == NULL) 
fatal ("alloc Term(): out of memory\n"); 
pi(n-1]).coeff - coeff; 
p(n-1l}.B = alloc_bound(); 
) 


return (p); 


Bound *alloc_ bound () 


/* 


———S—eSe@3w G3 ewe eS ee eee oe oe oe oe oe em oe oe wt wt oe oe om om et wt wt wt wt wt wt we im wm wm wm wt wt ee et ee es eee ee ee ee ee ee 


-— oe om oe ow om om ow et et et es ee oe om om om © om om om ot om ae et ee et ee ee we oe om wt oo oe oe oe oe we we ee = oe oe = Oe Oe ge Oe ee oo oe Oe oe ee oe = oe 


tfunction: 


- Allocate space for EF orig.nvar bounds entries and initialize 


each bound to -1,f orig.radix-1. 

- If p is NULL, allocate new space 
:globals: 

£ orig 
treturns: 

- A pointer to the Bound array 


Bound *p; 
char *malloc({); 
register i; 
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ee ee ee ee ee eee ee ee ete ee ee ee ee ee ee ee ee ee oe oe oe oe oe ee ee oe 


Lf ( (p= (Bound *)malloc(sizeof (Bound) *(F_orig.nvar))) == NULL) 
fatal ("alloc bound (): out of memory\n"); 


for (i-0; i < E orig.nvar; i++) [ 
pli]J.lower = -1; 
p(iJ.upper = E_orig.radix-1; 
) 


return(p); 
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init expr () 


:function: 

- Initialize E_ work, F_orig and F_ final 
7Jide Effects: 

E work 

E orig 

E final 


cc ee ee ee ce ce ce ge ee ee ce ee ee ee ee ee me ee ee ee ee ee ee ee ee ee we ee ee ee ee oe em we oe Oe oe oe om Om oe mm oe om oe oe 


EF work.T = NULL; 

F orig.T = NULL; 

EF orig.nvar = 0; 

F orig.nterm = 0; 

FE orig.radix = 0; 

FE final(0].T = NULL; 
ye final(1).T = NULL; 


dealloc expr (e) 
Expression ‘*%e; 


:function: 
- Deallocate the expression pointed to by e 


ee eee ee ee ee ee ee ee ee ee ee ee ee we me ee ee ee em ee ee ee ee ee ee em me ee we = 


term yi as 
register i? 


if (e->T f= NULL) [{ 
for (p = e->T,i-0; i < e->nterm;: i++) 
if (p{i].B t= NULL) 
free (p{i].B); 
free (p); 
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dup expr(E dest,E src) 
Expression hides. fuscc, 


:function: 
- Duplicate the expression pointed to by E sre by allocating 
and copying into the expression pointed to by E_dest. 
:calls: 
alloc _bound() 


wee ee een mcm nee ea eens eee eee on eea sea] =] ee a/ 
( 

Term *Te 

Bound tal ae 

register 2 a |e 

char *malloc(); 


E dest->radix = E src->radix; 
E dest->nvar e FE sre->nvar; 
E_dest->nterm = ESce--neerm, 


if ((T=(Term *)malloc (sizeo£ (Term) *(E dest=>nterm) )) == 00) 
fatal ("dup _expr(}: out of memory\n") ; 


for (1-0; i < Eusre-snterm: i117) | 
Tl{i}.coeff = £ sre->Tli}.coeff; 
T(i) 8 = alloc tbound() ; 
for (j-0; j <E srce->nvar;: jt) | 
Tli}.B[j}) .lower =~ E sre->T({i}.B(j) lower; 
T{i}).B(j).upper = E_sre->T(i).B(j). upper; 
} 


) 
E dest<>T = T; 
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/* mvitest.c 


- A test source file generator 


a / 


finclude <stdio.h> 

char 

‘version = "\ntiVLTEST V1.0 Copyright 1966 Naval Postgraduate School\n\n", 
*‘usage = “usage: mviltest {-sN -rN ~iN -tN -oN}\n", 


*help(} = { 
" -si! ~ Cycle the random seed N times (default = 0)", 
S -rH ~ Set radix to N (default = 4)", 
" -inN - Set number of inputs to N (default = 2)", 
° -tN ~ Set number of terms to N (default = 1)"*, 
is ~oN ~ Set number of output to N (default — 1)"*, 
NULL 


)? 


main(argce, argv) 
char *argv(}; 
( 
int nterms = 1; 
int reps = 1, 
radix = 4, 
nvar = 2, 
éoceff = 1; 
int seed = 0; 
int 4,4,k; 
int *b, *get hound (); 
float random(); 
char ol airs 
register ap; 


for (ap<1; ap < argc; aptt) [ 
p = argv[ap); 
if (*ptt — ’=") f{ 
while (*p) { 
switch (ptt) 
case ‘r’': 
sscanf (p, "$d", Gradix); 
Spm ONO’ s 
break; 
case ‘i’: 
sscanf (p, "$d", &nvar) ; 
Deze NO"? 
break; 
case ‘t's: 
sscanf (p, "$d", &nterms) ; 


‘p= *\0'; 
break; 
case 'o’: 
sscanf (p, "td", &Greps); 
‘pose \0Oes 
break; 
case 's': 
sscant (pop, ta ,;6s¢eq)), 
*p & "\0';: 
break; 
case '-': 
fprintf(stderr, "%$s%s", version, usage) ; 
for (ap=0; help(ap); apit) 
fprintf(stderr, "&s\n",help[ap)); 
exit (0); 
default: 
fprintf(stderr, "mviltest: unknown switch - %c\n",p[(-1]}); 
fprintf(stderr, "%3s", usage); 
exit(1); 


J L} 


) 
for (1-0; 4 < seed; i++) random(}; 


for (1 = 0; 1 < reps; itt) [ 
printf ("%td:%d:", radix,nvar); 
for (j=0: } < nmterms; 4jt+) ( 
printf ("\n 4%da", (int) (random() * (float) (radix-1))+41); 
for (k=1: k <= nvar; kt) ( 
be get bound (radix); 
printf ("*x%d(%d, $d) ",k,b(0),b(1)); 
) 


) 
printf (";\n"); 


int *get_ bound (radix) 


( 
static int b[3]; 
float random (); 


for (77) { 
b(O}) = (int) (random() * ((float) (radix)-0.001)); 
b(1) = (int) (random() * ( (float) (radix)-0.001)); 
if (b(0) <= b{1)) 
break; 
) 
return(b); 


) 


float random () 


( 
static long a=100001; 


a= (a*125) %& 2796203; 
return ((float)a/2796203); 
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f Makefile for mvl 

t 

CF LAGS=-O 

OBJS=main.o dm.o pa.o common.o parse.o alloc.o 


all: mvl mvitest 


mvl: $(OBJS) 
cc $(CFLAGS) $(OBJS) -o mv1l -1Ll 


mvitest: mvitest.o 
c&é $(CFLAGS) mvitest.o -o mvitest 


main.o: main.c | 
Aalloc.o: alloc.c 
parse.o! parse.c 
dm.o: dm.c 
pa.o!: pa.c 
common.o: common.c 
parse.c: parse.1 
lex parse.l 
mv -f lex.yy.c parse.c 


mvitest.o: mvitest.c 
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